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This  work  focuses  on  two  fundamental  problems  in 
scattering  theory.  The  first  is  centered  around  the 
electronic  basis  used  to  expand  the  solution  of  the  time 
independent  Schrbdinger  Equation.  The  other  consists  of 
how  to  handle  the  problem  for  the  nuclear  degrees  of 
freedom  when  more  than  one  electronic  basis  function  is 
included  in  the  expansion. 

In  chapter  two  it  is  shown  that,  in  the  adiabatic 
electronic  representation,  the  Hamiltonian  is  not 
hermitean  if  more  than  one  electronic  basis  function  is 
included  in  the  expansion.  It  is  painted  out  that  in  such 


cases  alternative  electronic  representations  may  be  more 
suitable.  The  adiabatic  and  diabatic  representations  are 
reviewed  and  the  "nearly  adiabatic"  representation  and  a 
representation  based  on  a mini/tiieation  procedure  are  intro- 
duced. Calculations  are  done  to  compare  the  results  from 
the  diabatic  representation  and*  the  representation  arising 
from  the  minimization  procedure. 


Chapter  three  considers  the  problem  for  the  nuclear 
degrees  of  freedom  in  multi-surface  systems.  A "common" 
eikonal  is  used  for  the  nuclear  wave  function  satisfying 
the  time  independent  Schrfidinger  equation.  This  "common" 
eibonal  is  obtained  from  a modified  Hamllton-Jacobi  equation 
that  involves  an  average  potential.  Implementing  the  short 
wavelength  approximation  and  transforming  to  a time 
dependent  picture  leads  to  a set  of  first  order  differential 
equations  in  time  that  determine  the  expansion  coefficients. 
It  is  shown  that  these  equations  along  with  the  differential 
equations  obtained  from  the  Hamilton-dacobi  like  equation 

first  order  differential  equations  in  time  thet  are  formally 
equivalent  to  Hamilton's  equations  of  motion. 

This  procedure  is  applied  to  the  collinear  H3'"  system. 
Hyparspherical  coordinates  are  used  for  the  nuclear  degrees 
of  freedom  and  the  adiabatic  potential  energy  surfaces  and 


electronic  coupling  terms  are  obtained  from  the  method  of 
Diatomics  in  Molecules.  A transformation  to  the  diabatic 
representation  1s  made  and  the  trajectories  are  celculated 


viii 


in  thiB  representation.  Typical  trajectories  are  presented 
in  chapter  six  and  a comparison  o£  the  total  electron 
transfer  probability  with  those  from  a quantum  mechanical 
study  is  made.  The  results  are  encouraginq. 


CHAPTER  1 
INTRODUCTION 


1-1  General  Problaw 

This  study  priirarily  addresses  the  problem  of  molecular 
reactions  where  more  than  one  electronic  state  is  energeti- 
cally accessible.  These  processes  are  usually  referred  to 
as  non-adiabatic  collisions.  The  collisions  are  assumed  to 
take  place  between  combinations  of  atomic  or  molecular 
fragments  that  can  be  either  charged  or  neutral.  Since  the 
origins  of  the  formalism  to  be  presented  are  based  in 
quantum  mechanics,  a few  comments  on  the  general  quantum 
mechanical  treatment  of  the  problem  seem  appropriate. 

The  general  solution  to  the  Schrbdinger  equation 
involving  both  nuclear  and  electronic  degrees  of  freedom  is 
a function  of  the  nuclear  and  electronic  coordinates.  Until 
programs  are  available  to  numerically  solve  these  multi- 
dimensional equations,  some  form  of  approximation  to  the 
solution  must  be  made.  Generally,  the  solution  is  expanded 
in  an  electronic  basis.  This  in  itself  does  not  simplify 

the  electronic  problem  can  be  solved  independently  of  the 
dynamics  of  the  nuclei. 

The  approximation  referred  to  above  is  the  well-hnown 
Born-Oppenheimer  approximation  (Born  and 


Oppenheimer,  19271. 


approximation  is  physically  based 


electron's  mass  is  orders  of  magnitude  smaller  than  the  mass 
of  a nucleus  and  hence  that  the  electronic  motions  are  much 
faster  than  nuclear  motions.  Then  one  would  expect  that  the 
dependence  of  the  wave  function  describing  electronic  motion 
on  nuclear  variables  could  be  neglected.  In  other  words, 
this  approximation  assumes  that  the  nuclear  kinetic  energy 
operator  can  treat  the  electronic  wave  functions  as  con- 
stants and  leads  to  solutions  of  the  electronic  problem  that 
depend  only  on  the  positions  of  the  nuclei. 

with  this  approximation,  the  electronic  problem  is 
solved  with  the  positions  of  the  nuclei  held  fixed.  Since 
the  electronic  problem  mast  be  solved  for  each  internuclear 
configuration,  the  electronic  functions,  energies,  etc.  are 
said  to  depend  parametrically  on  the  nuclear  coordinates. 
Since  the  dependence  of  the  nuclear  wave  function  on  the 
electronic  coordinates  has  been  ignored,  the  solution  to  the 
problem  is  in  the  form  of  a linear  con^lnation  of  products 
of  electronic  functions  that  depend  parametrically  on  the 
nuclear  coordinates  and  nuclear  functions  that  depend  only 
on  the  nuclear  variables. 

This  approximation  is  an  extreme  simplification  to  the 
most  general  solution  in  that  the  electronic  problem  is 
solved  without  a full  dynamical  knowledge  of  the  nuclei.  In 
other  words,  this  approximation  assumes  that  the  electronic 
problem  can  be  solved  independent  of  the  nuclear  velocities. 
This  approximation  may  not  be  valid  when  the  nuclear  speeds 


are  high  (comparable  to  those  of  the  electrons) . But  in 
this  study,  the  nuclear  speeds  will  be  slow  compared  to 
electronic  speeds  so  that  this  approximation  will  be 
assumed  to  hold  to  a high  degree  of  accuracy. 

From  the  preceding  discussion,  it  is  apparent  that  there 
are  two  essential  steps  to  salving  the  problem.  The  first 
is  to  choose  an  electronic  basis  that  is  convenient  for  the 
calculation  and  to  solve  the  electronic  problem  within  the 
Born-Oppenheimsr  approximation.  This  will  be  the  subject 
matter  of  the  next  section  and  chapter  two.  The  second  step 
is  to  expand  the  solution  in  this  basis  and  solve  for  the 
nuclear  expansion  coefficients.  This  will  be  the  topic  of 
the  last  section  of  this  chapter  and  chapter  three. 

1~2  Electronic  Hotions 

is  was  discussed  in  the  previous  section,  the  solution 
of  the  full  Schrodinger  equation  is  expanded  in  an  electronic 
basis  and  the  electronic  problem  is  solved  with  the 
nuclei  held  fixed.  In  general,  the  expansion  in  a set  of 
n electronic  states  will  result  in  a square  nxn  matrix 
repreaestation  of  the  electronic  Hamiltonian,  i.e.  the  full 
Hamiltonian  without  the  nuclear  kinetic  energy  operator. 
Couplings  between  the  nuclear  expansion  coefficients  will 
occur  through  the  off-diagonal  matrix  elements  of  the  elec- 
tronic Hamiltonian  and  through  the  electronic  matrix 
representation  of  the  nuclear  kinetic  energy  operator. 


all  of  these  couplings  and 


brief  discussion  of  some  of 


the  more  familiar  representatione  and  their  properties 

One  of  the  moat  videly  used  electronic  representations 
is  the  adiabatic  representation  (Bom  and  Oppenheimec,  1927; 
Born  and  Huang,  19S4I.  This  representation  has  both  been 
extensively  used  in  bound  state  calculations,  i.e.  in  calcu- 
lations where  the  dynamics  of  the  nuclei  can  be  ignored,  and 
in  scattering  calculations.  It  is  essentially  characterized 
by  giving  a diagonal  matrix  representation  of  the  electronic 
Hamiltonian. 

The  (diagonal)  matrix  elements,  often  referred  to  as 
eigenvalues,  of  this  representation  give  rise  to  potential 
energy  hypersurfaces.  For  adiabatic  reactions,  the  reac- 
tants and  products  are  both  on  tihe  sane  potential  energy 
hypersurface  while  for  non-adiabatic  reactions  the  surface 
for  the  reactants  is  different  from  that  of  the  products,  h 
conseguence  of  this  representation  is  that  surfaces  of  the 
same  electronic  symmetry  do  not  cross  (Moiseiwitsoh,  19«1), 
i.e.  oouplings,  in  most  cases,  do  not  occur  through  the 
electronic  Hamiltonian. 

These  surfaces  are  essential  in  many  formulations  of 
the  problem.  As  mentioned  in  the  previous  section,  the  eigen- 
values of  this  representation  depend  parametrically  on  the 
nuclear  coordinates  so  that  the  electronic  problem  must  be 
solved  for  each  nuclear  configuration.  Since  scattering  cal- 
culations normally  require  many  nuclear  configurations,  an  ab 
initio  treatment  of  the  electronxc  problem  would  almost  surely 
be  prohibitive.  This  offers  motivation  and  support  for  using 


more  approximate  procedures  auoh  as  Diatomics  in  Molecules 
(Ellison,  1963)  or  Polyatomiea  in  Moleculea  (Tully,  1977) 
which  mahe  uae  of  either  experimental  or  theoretical  know- 
ledge of  the  diatomics  or  polyatomics  to  construct  the 


surfaces. 

Since  the  electronic  Hamiltonian  is  diagonal  in  this 
representation,  the  couplings  between  the  nuclear  expansion 
coefficients  arise  from  the  matrix  representation  of  the 
nuclear  )tinetlo  energy  operator  and  gradient  in  this  basis. 

The  inclusion  of  these  terms  in  a scattering  calculation  leads 
to  two  main  difficulties,  first,  efficient  computer  programs 
are  not  available  for  solving  sets  of  equations  that  contain 
these  terms  and  second,  as  has  been  noted  (Smith,  1969),  it 
all  coupling  terms  are  included,  the  Hamiltonian  in  this 
electronic  basis  is  not  herraitean.  Although  the  first  dif- 
ficulty is  a computational  one  that  could  be  conceivably 
overcome, the  second  difficulty  introduces  theoretical  problems 
such  as  nonconservation  of  energy  and  flux  or  the  need  to 
introduce  a second  set  of  nuclear  coefficients  that  satisfy 
the  adjoint  operator  of  the  original  problem. 


is  well  understood 
tations  don't  have 
related 


:an  be  surmised, this  representation  is  probably  not 
useful  for  processes  that  exhibit  non-adiabatlc 
This  representation  is  however  still  essential 
he  solution  of  the  electronic  problem  in  this  basis 
theoretically  sound.  Other  represen- 
lave  this  property  so  that  they  must  be 
adiabatic  basis  by  way  of  unitary 


transComations  in  order  to  fom  a matrix  representation  of 
the  electronic  Uamiltonian.  That  is,  the  eigenvalues  of 
the  adiabatic  representation  are  convenient  to  form  a matrix 
representation  of  the  electronic  Hamiltonian  in  any 
electronic  representation. 

There  are  a variety  of  representations  other  than  the 
adiabatic  one  that  can  be  used  in  a scattering  calculation. 
They  are  based  on  either  completely  or  partially  eliminating 
the  couplings  due  to  the  nuclear  kinetic  energy  operator  or 
gradient  or  on  minimising  the  couplings  in  a way  to  be 
discussed  later.  These  schemes  have  one  thing  in  common  in 
that  coupling  terms  present  in  the  adiabatic  representation 
that  are  eliminated  or  minimised  are  replaced  by  electro- 
static terms  that  give  rise  to  a non-diagonal  matrix 
representation  of  the  electronic  Hamiltonian. 

The  electronic  problem  will  be  the  subject  matter  of 
chapter  two.  Four  representations  including  the  adiabatic 
one  are  discussed  and  some  of  their  advantages  and  disadvan- 
tages are  pointed  out.  It  will  be  shown  that  care  must  be 
used  in  selecting  the  representation  that  would  be  most 
suitable  for  a scattering  calculation. 

1-3  nuclear  Motions 

Once  the  electronic  representation  has  been  chosen 
there  remains  the  problem  of  solving  for  the  nuclear  expan- 
sion coefficients.  Generally  they  are  solved  for  quantum 
mechanically,  semi-clasalcally  or  classically.  Since  quan- 
tum mechanical  and  semi -classical  treatments  of  this 


discussion 


problem  are  not  topics  of  this  study,  only  a brief 
of  these  methods  follows.  Comments  on  these  subjects  will 
be  restricted  to  be  general  and  provide  i^tivation  for 
using  the  method  to  be  developed. 

Since  the  electronic  problem  was  treated  quantum 
mechanically,  it  would  seem  natural  to  treat  the  nuclear 
problem  quantum  mechanically  too.  Indeed,  this  does  give 
dynamically  consistent  formalisms  but  all  quantum  mechanical 
approaches  to  this  problem  have  a serious  disadvantage. 

Since  the  nuclei  have  vibrational,  rotational  and  trans- 
lational degrees  of  freedom,  a quantum  mechanical  approach 
even  for  inelastic  collisions  would  require  expanding  the 
nuclear  expansion  coefficients  in  at  least  internal  states 
corresponding  to  the  vibrational  and  rotational  degrees  of 
freedom.  Since  normally  many  internal  states  are  energeti- 
cally accessible,  this  expansion  leads  to  a large  number  of 
coupled  differential  equations  to  be  solved  simultaneously. 
This  is  actually  one  of  the  simpler  cases  In  that  if  a 
partial  wave  expansion  is  made  a normally  much  larger  set  of 
equations  is  obtained. 

The  problem  becomes  even  more  difficult  for  reactive 
collisions.  In  this  case  there  is  more  than  one  asymptotic 
Hamiltonian  to  be  considered.  Expansions  in  internal  states 
for  each  asymptotic  Hamiltonian  must  be  made.  Not  only  does 
the  number  of  coupled  equations  increase  but  basis  elements 
corresponding  to  different  asymptotic  Hamiltonians  are  not 


necessarily  orthogonal  and  can  lead  to  problems  such  as 
overcoirpleteness. 

It  would  seem  then  that  even  though  a quantum 
mechanical  approach  is  theoretically  appealing,  it  leads  to 
some  very  fundamental  difficulties.  Even  if  approximations 
were  introduced  to  reduce  the  number  of  coupled  equations, 
it  would  be  helpful  to  introduce  a formalism  that  would 
avoid  these  difficulties  as  much  as  possible.  Before 
pursuing  this  however  a few  comments  on  the  semiclassical 
approach  to  non-adiabstic  collisions  without  nuclear 
rearrangements  will  be  made. 

Semi -classical  approaches  to  non-adiabatic  processes 
without  nuclear  rearrangemeht  have  been  used  extensively  and 
many  good  reviews  are  available  on  this  subject  (KcDowell 
and  Coleman,  19701.  These  approaches  normally  start  with 
the  time  dependent  Schrodinger  equation  for  the  electronic 
Hamiltonian  and  assume  that  the  trajectories  of  the  nuclei 
are  known.  One  then  needs  to  solve  a set  of  first  order 
differential  equations  in  time  for  the  expansion  coeffi- 
cients of  the  electronic  wave  function.  Simultaneous 
integration  of  these  equations  gives  these  coefficients  as  a 
function  of  time  and  their  values  at  the  final  tine  are 
related  to  the  probability  of  non-adiabatic  transitions  for 
the  collieion. 

As  mentioned  in  the  preceding  paragraph,  the  nuclear 
trajectories  in  these  approaches  are  assumed  known.  These 


different  approximations  such  as  the  Impact  parameter 


method  (hi9h  velocity  approxinationl  or  the  perturbed 
stationary  state  method  (low  velocity  limit).  Besides  the 
limitations  caused  by  not  determining  the  trajectories  from 
ab  initio  considerations  the  extension  of  these  approaches 
to  include  nuclear  rearrangement  is  not  trivial. 

Since  the  masses  of  the  nuclei  are  so  much  heavier  than 
those  of  the  electrons,  it  would  seem  possible  that  some  or 
all  of  the  nuclear  degrees  of  freedom  could  be  treated 
classically.  This  approach  has  been  successful  for  the  case 
of  one  potential  energy  surface,  i.e.  in  the  absence  of 
couplings  between  nuclear  and  electronic  degrees  of  freedom. 
For  such  systems,  the  nuclei  evolve  according  to  classical 
equations  of  motion  on  a quantum  mechanical  potential  energy 
surface.  The  initial  conditions  of  the  trajectories  are 
obtained  by  talcing  a Monte  Carlo  sampling  of  the  possible 
initial  states  of  the  system.  This  approach  is  however  not 
as  straightforward  when  more  than  one  surface  is  included. 

A primary  concern  in  this  study  is  to  develop  a 
formalism  that  extends  the  classical  treatment  of  nuclear 
degrees  of  freedom  to  systems  that  are  characterized  by  more 
than  one  electronic  state.  This  avoids  the  difficulty  of 
expansions  in  internal  states  required  by  quantum  mechanical 
treatments.  Also  assumptions  about  the  nuclear  trajectories 
are  not  needed  since  they  are  determined  from  classical 
equations  of  motion.  This  also  has  the  advantage  that 
individual  trajectories  give  a clear  conceptual  picture  of 
the  collision  event. 


Briefly,  the  approach  to  be  used  I 
independent  Schrodin9er  equation. 


wave  function  ii 


expanded  in  a set  of  electronic  states.  The  nuclear  expan- 
sion coefficients  are  written  as  products  of  an  amplitude 
and  a common  phase  which  is  proportional  to  the  eikonal. 

The  gradient  of  the  eikonal  is  required  to  satisfy  the 
Bamilton-Jacobi  equation  whose  potential  is  determined  from 
the  quantum  mechanical  equations.  Upon  implementing  the 
short  wavelength  approximation  and  making  a transformation 
into  time, it  is  found  that  the  nuclear  expansion  coeffi- 
cients satisfy  first  order  differential  equations  in  time 
and  that  the  gradients  of  the  eikonal  become  the  classical 
momenta  of  the  nuclei.  If  a convenient  form  for  the  nuclear 
expansion  coefficients  is  chosen,  it  is  found  that  not  only 
the  nuclear  positions  and  momenta  but  also  the  expansion 
coefficients  satisfy  Hamilton^  equations  of  notion. 

Although  other  approaches  (Meyer  and  Hiller,  19791  have 
obtained  similar  results, the  author  considers  this  treatment 
to  be  on  a more  sound  theoretical  foundation.  The  develop- 
ment of  thie  forinellsm  will  be  the  subject  matter  of 
chapter  three. 

The  formalism  of  chapters  two  and  three  is  applied  to 
the  collinear  H}  system.  This  system  was  chosen  partly 
because  of  the  presence  of  large  non-adiabatic  coupling 
terms.  Also,  due  to  the  relatively  light  masses  of  the 
nuclei,  this  system  should  provide  a good  test  of  the 
theory.  It  is  also  interesting  from  the  viewpoint  that  not 


inelastic  processes 


only  do  elastic  and 
exhibits  reactions  and  or  rearrangements.  Finally  is 
adequately  described  with  two  electronic  states  which  Is  a 
natural  starting  point  for  an  application  of  the  theory. 

The  electronic  and  nuclear  parts  of  the  problem  for  Hj 
will  be  developed  in  chapters  four  and  five  respectively. 

The  electronic  problem  will  be  solved  with  the  method  of 
Diatomics  in  Molecule  and  the  nuclear  problem  will  be 
solved  in  hyperspherieal  coordinates.  This  convenient 
choice  of  coordinate  system  will  be  discussed  at  length 
in  Appendix  one. 

This  study  will  conclude  with  a presentation  of  the 
calculations  in  chapter  six.  A brief  background  on  previous 
experimental  and  theoretical  studies  for  the  Uj''  system  will 
also  be  given.  A comparison  with  a quantum  iKchanical 
calculation  (Top  and  Baer,  1977)  for  collinear  Hj*  will  be 
made  and  this  work  will  close  with  some  comments  and 
conclusions . 


CHAPTER  2 

ELECTRONIC  REPRESENTATIONS 

2-1  Introduction 

The  adiabatic  approximation  has  played  a central  role 
in  the  study  of  molecular  processes  since  its  introduction 
by  Born  and  Oppenheimer  (Born  and  Oppenheisiec.  1927;  Born 
and  Huang,  1954).  The  approximation  introduces  a basis  of 
electronic  states  that  provide  an  adiabatic  representation 
for  electronic  operators.  Non-adiabatic  collisions, 
originally  studied  by  Landau,  Zener  and  stueckelberg, 
(Landau,  1932;  Zener,  1932 ; Stuec)celberg , 1932:  Nikitin, 

1970)  require  information  on  the  momentum  couplings  of  the 
adiabatic  representation.  Many  models  have  been  developed 
to  incorporate  the  couplings  (child,  1979;  Tully,  1976: 
Garrett  and  Truhlar,  1980;  Delos,  1981). 

In  cases  where  the  adiabatic  couplings  can  not  be 
ignored,  a different  electronic  basis  and  corresponding 
representation  itiay  prove  to  be  useful.  One  such  alternative 
is  a diabatic  representation  (Smith,  1969)  defined  so  that 
momentum  couplings  are  exactly  eliminated  and  transitions 
occur  only  through  Coulomb  interactions  of  electrons  and 
nuclei.  This  is  done  by  introducing  the  eigenstates  of 
the  momentum  operator.  In  its  original  version  this 
representation  was  oriticited  because  the  electronic  states 


of  the  new  representation  could  not  change  with 
interiDolecular  distance,  except  for  phase  factors  (Gabriel 
and  Taulbjerg,  1974).  Working  however  with  finite  bases  one 
can  define  a diabetic  representation  by  requiring  that  the 
matrix  of  the  moinentuin  operator  is  zero.  This  introduces 
a matrix  unitary  transformation  from  the  adiabatic  to  the 
diabetic  set  of  states  which  does  change  with  intermolecular 
distance.  Given  this  matrix,  one  can  transform  the  matrix 
of  Coulomb  interactions  to  the  new  basis. 

Although  this  procedure  is  mathematically  rigorous,  it 
may  lead  to  complications  in  the  physical  description  of 
collisional  processes.  As  we  shall  see,  depending  on  the 
magnitude  of  the  momentum  couplings  in  the  origihal  adia- 
batic representation,  the  diagonal  eles^nts  of  the  Coulomb 
interaction  matrix  in  the  new  diabatic  representation  may  be 
far  from  physically  meaningful.  For  large  iiiomentum  coup- 
lings these  diagonal  elements  may  repeatedly  cross ; while 
for  small  momentum  couplings  they  may  be  far  from  adiabatic 
potentials  in  regions  where  these  are  physically  meaningful. 
The  latter  problem  can  be  particularly  significant  in 
studies  of  reactive  atom-diatom  collisions  because  the  new 
representation  may  give  the  correct  description  of  the 
reactant  potentials  but  a completely  unphysical  one  for  the 
product  potentials. 

These  difficulties  result  from  using  finite  bases  and 
from  the  differential  equation  satisfied  by  the  unitary 
transformation,  which  is  of  first  order  in  the  intermolecular 
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position  variables.  It  follows  that  given  the  knovm 
boundary  conditions  for  large  distances  the  transformation 
and  diabatic  potentials  are  mathematically  determined  for 
all  shorter  distances,  which  does  not  leave  any  room  for 
physical  considerations. 

The  aims  of  the  present  chapter  are  to  introduce  a new 
diabatic  representation  which  leads  to  physically  well 
behaved  potentials  <in  the  sense  to  be  described),  and  to 
show  how  it  is  constructed  around  pseudocrossings.  It 
starts  with  the  adiabatic  potentials  and  couplings  which 
are  obtained  in  eleotronic  structure  calculations.  The 
procedure  is  based  on  the  minimisation  of  coupling  terms,  and 
provides  a criterion  to  determine  the  range  of  kinetic 
energies  over  which  it  is  justifiable  to  neglect  couplings 
altogether. 

Other  diabatic  representations  are  possible  and  have 
been  introduced  by  means  of  physical  arguments  (o'Halley, 
1967).  For  energetic  atom-atom  collisions,  several  of  the 
representations  have  been  extended  to  incorporate  electron 
translation  factors  in  order  to  satisfy  asymptotic  condi- 
tions (Delos  and  Thorson,  1979;  Delos,  19B1).  These 
extensions  shall  not  be  considered  because  the  immediate 
aims  refer  to  therroal  and  hyperthermal  collisions.  Adia- 
batic and  diabatic  representations  have  also  been  introduced 
for  atom-diatom  collisions  (Baer,  1975!  Top  and  Baer,  1977a). 

Numerical  studies  of  electronic  states  in  various 
representations  include  calculations  of  potentials  and 
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their  couplings  for  atom-atom  (Redmon  and  Hicha,  1974; 
Kimrich  and  Truhlar,  1975;  Evans,  Cohan  and  Lane,  1971) 
and  atom-diatom  systems  IRebentrost  and  Lester,  1977; 

Tally,  19801.  A great  deal  of  related  worJc  has  also  been 
done  on  diabetic  molecular  orbitals  and  their  energies 
ILichten,  1963;  Briggs,  1976;  HcCarroll,  1976),  to  which  the 
developments  in  this  chapter  could  also  be  applied. 

To  illustrate  sotoe  of  the  numerioal  aspects,  some 
results  are  briefly  mentioned  for  Bj  and  R));  in  the 
collinear  conforsiations.  For  a basis  of  two  electronic 
states  (the  two  lowest  states  of  symmetry  for  H;  ; the 
lowest  and  states  for  PH^),  the  transformation  from 
the  adiabatic  to  a diabatic  basis  depends  on  the  integral 


where  is  the  momentuio  coupling  in  the  nuolear  variables 
X,  and  D is  their  domain.  As  shall  be  seen  in  Section  2-4, 
the  standard  diabetic  representation  (Smith,  1969)  »or)ts 
well  when  the  integral  equals  tt/2  but  not  when  it  differs 
appreciably  from  7r/2.  For  atom-diatom  mass-weighted  Jacobi 
variables  (Z,a)  for  the  intemolecular  and  internal  coordi- 
nates, a transformation  to  polar  coordinates  r ■ (2®  + a*)** 
and  6 * tan  hzfl\  leads  to  values  of 


which  for  go  from  0.3  to  1.6  as  r varies  from  4 a.u.  to 
12  a.u.  (Tully,  1976) , while  for  FHj  they  stay  around  0.7 
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as  r varies  from  B a.u.  to  12  s.u.  (Tally,  19S0).  Other 
numerical  examplee  can  be  found  in  the  recent  literature 
on  adiabatic-diabatic  transforroations  for  slow  nuclear 
degrees  of  freedom,  where  the  Integrals  are  instead  larger 
than  n/2  and  multiple  crossings  occur  (Baer,  Drolshagen  and 
Toennies,  1980). 

Given  the  wide  variety  of  problems  where  pseudo* 
crossings  may  occur,  this  chapter  shall  not  concentrate  on 
a given  physical  system  but  shall  instead  construct  a model 
of  potentials  and  couplings  with  parameters  that  will  be 
varied  around  physical  values.  The  shape  of  the  potentials 
and  couplings  are  similar  to  those  calculated  (Tully,  197B) 
for  H:  and  the  physical  parameters  relate  also  to  this 

system. 

The  adiabatic,  diabatic  and  nearly  adiabatic  represen- 
tations will  be  briefly  developed  in  Section  2-2.  This  will 
be  followed  in  Section  2-3  with  a general  development  of  the 
minimization  procedure.  Section  2-4  will  give  a detailed 
treatment  of  the  two  electronic  state  problem  in  one 
dimension.  A comparison  of  the  results  for  two  of  the 
diabatic  procedures  will  be  given  in  Section  2-5  and  the 
chapter  will  close  with  a discussion. 

2-2  Electronic  state  Representations 

He  consider  to  begin  with  a molecular  system  with  n 
nuclei,  in  a center-of-mass  reference  frame.  Introducing 
cartesian  coordinates  and  the  nuclear  position  vectors 


i*l  to  n},  the  nuclear  kinetic  energy  operator  ie 
expreased  aa 

^nu  * ■ ^ 3’/3yJ  (2-1) 

where  nt^  is  the  mass  of  the  1^  nucleus  and  we  have  used 
atomic  units  (fi*l) , Introducing  a change  of  variahles, 


(2-2) 


total  masS(  gives 


Defining  next  a nuclear  m 
with  n orthogonal  components  and 


(2-3) 


(2-4) 


where  X = (x^,  , . . x^),  the  kinetic  energy  operator  (becomes 

^nu  * (2M)'l  Pjjjj  . . -(2M)*l3V3X*  . (2-S) 

For  the  special  case  of  ayatema  with  two  nuclei,  the 
roomentun  operator  in  the  center-o£-masa  coordinate  system 


(2-6) 


where  R la  the  relative  position  vector  between  the  nuclei. 
For  the  three-nuclei  ayatein,  mass  weighted  Jacobi  coordi- 
nates also  satisfy  the  above  conventions. 

Solutions  of  the  time  independent  Schrddinger  equation 

(H-E)|’MX)>  . 0 , (2-7) 


H - T^ylS/3X)  ♦ H^^(X)  . (2-8) 

and  the  nuclear  positions  X were  defined  above  and 
is  the  electronic  Hamiltonian  including  nuclear-nuclear 
repulsion  terms.  The  brac)iet  notation  refers  to  the 
electronic  coordinates  and  involves  an  integration  only  over 
electronic  coordinates.  Invo)(ing  a separation  of  electronic 
and  nuclear  variables  leads  to  the  solution  having  the 


|t(X)>  - Z|». (X)>»^(X)  (2-9) 

where  one  requires  that  the  electronic  functions,  form 
a complete,  orthononDal  set  at  each  X. 


SubBtituting  Eq.  (2-9) 


Eq.  (2-7)  and  multiplying  Ccom  the 


left  by  <9^ (X) | leads  to 


(2-10) 


where  operates  on  all  factors  to  the  right.  The 
development  thus  far  is  completely  general  but  not  of  much 
use.  Hore  useful  representations  may  be  obtained  by 
requiring  the  electronic  functions  to  obey  additional 
properties  besides  those  of  completeness  and  ortho- 
normality and  these  will  be  briefly  discussed  below. 

2-2a  The  Adiabatic  Representation 

This  highly  useful  and  widely  used  representation 
(Bom  and  Oppenheimer,  1927;  Born  and  Huang,  19S4)  is  based 
on  the  requirement  that  the  electronic  functions  satisfy 


where  the  superscript  "a"  denotes  the  adiabatic  representa- 
tion. Using  this  basis  in  Eq.  (2-10)  leads  to 


where  operates  on  all  factors  to  the  right.  As  is  well 
known,  adiabatic  potential  surfaces  v5.  of  the  same  symmetry 


(X) 


(2-11) 


S<«®|T^„|«®>*®(X) 


(2-12) 


all  i,j 


n point  would  be  This  basis  can  b 


[|^(X„)>=  |*=<X)  A(X,X^, 

[♦''>  k=d,a  and  A is  a nxn  matrix. 

Eq.  (2-13}  one  finds  that  the  transformation  must 
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.“lo  - ™«i.. 


(2-11). 


approximation  would  be  valid  in  these  regions.  In  region 
lb)  the  potentials  do  approach  each  other  relatively  closely 
so  the  adiabatic  approximation  would  no  longer  be  valid  and 
one  should  either  include  the  momentum  coupling  terms  in 
the  adiabatic  basis  or  use  a diabatic  representation.  Then 
physically  well  behaved  diabatic  potentials  should  fulfill 
the  following  conditions.  In  region  (c),  as  shown  in  Fig. 
(1).  the  diagonal  elements  of  the  diabatic  potential  matrix 
vfi  and  Vji  should  coincide  with  the  adiabatic  potentials 
vfi  and  vfi,  respectively.  In  region  (a)  vfi  and  vti  should 
coincide  with  vfj  and  v?i,  respectively.  In  region  (b)  the 
diagonal  elements  of  the  diabatic  potential  matrix  should 
vary  smoothly  and  exhibit  a single  crossing.  Finally,  the 
off  diagonal  matrix  elements  V12  and  vfi  should  vanish  in 
regions  (a)  and  (c)  far  from  the  crossing. 

As  is  seen  in  Eg.  (2-15),  the  transformation  matrix 
A(X,X^)  satisfies  a first  order  differential  equation  and 
one  therefore  has  only  one  boundary  condition  at  one’s 
disposal.  This  will  ensure  proper  behavior  in  one  of  the 
regions  in  Fig.  (1),  usually  chosen  to  be  region  (e),  but 
the  behavior  in  the  other  regions  will  depend  on  the 
coupling  terms  ***|Pnul**’  so  that  this  procedure  does  not 
in  general  guarantee  diabatic  states  that  satisfy  the  above 
conditions. 

2-2c  The  Nearly  Adiabatic  Representation 

As  was  pointed  out  in  Section  1-2,  one  of  the 
difficulties  that  arise  when  using  an  arbitrary  electronic 
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With  the  fom  of  given  by  &q.  (2-4)  one  seee  that  the 
tranaforroation  aatisfiee  a second  order  differential 
equation.  Since  there  are  two  boundary  conditions,  use  of 
this  representation  could  ensure  proper  behavior  of  the 
electronic  potentials  on  both  sides  of  the  pseudocrossinq. 

It  should  be  noted  that  since  second  order  partial 
differential  equations  are  usually  difficult  to  solve,  the 
usefulness  of  this  representation  may  be  somewhat  limited. 

It  should  however  provide  an  alternative  in  cases  where  the 
non-Hermitian  term  in  the  adiabatic  representation  can  not 
be  ignored  and  the  strictly  diabatic  representation  gives 
unphysical  results. 

2-3  The  Mininiaation  Procedure 

Physically  ill  behaved  diabatic  potentials  can  be  a 
serious  drawback  in  studies  of  scattering  processes.  The 
procedure  to  be  developed  here  starts  from  a different  point 
of  view.  Instead  of  completely  neglecting  the  momentum 
coupling  terms,  our  procedure  requires  well  behaved  diabatic 
potentials.  It  introduces  a unitary  transformation  with 
parameters  to  be  varlatlonally  chosen,  and  then  minimizes 
the  momentum  coupling  in  the  pseudocrossing  region. 

Since  the  electronic  Hamiltonian  is  usually  only  known 
in  the  adiabatic  representation,  a unitary  transformation 
between  the  representations  is  necessary.  In  this  procedure, 
a transformation  B will  be  chosen  to  guarantee  physically 
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so  that  Eg.  (2-261  beeomas 
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(2-34) 
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expression  should  be  helpful  In  determining  the  range  of 
the  relative  kinetic  energy  where  the  momentum  coupling 
can  be  neglected. 

Neglecting  the  x"  coupling  terms  should  also  be 
justified.  Defining 

^ (2-361 

and  assuming  the  adiabatic  basis  is  complete,  K*"  can  be 
related  to  K®  by  Eg. (2-28)  with  K®  obtained  by  letting 
®nu  P"  '**iPnul**^  give 

K®  ” (2-37) 

where  the  P®  • P*  term  cooes  from  the  completeness  of  *®. 
Since  the  literature  usually  only  reports  couplings  between 
the  lowest  electronic  states,  this  expression  may  not  be 
useful.  If  a valence  bond,  Diatoaics  in  Molecules  or  other 
method,  was  used  to  obtain  the  x-dependent  expansion 
coefficients  for  the  electronic  states,  the  K*  terms  could 
be  obtained  by  numerical  methods.  Once  K®  has  been  calcula- 
ted, it  is  straightforward  via  Eq.  (2-28) to  determine  k“. 

A justification  for  neglecting  it  would  depend  on  whether 


for  each 


examples  given  later  in  this  paper,  only  a functional  form 
for  P®  will  be  obtained  and  the  terras  k“  will  not  be 
considered.  Normally,  though,  if  the  coupling  from  is 
small,  one  would  expect  that  the  coupling  from  would 

2~<  "nie  Two  Electronic  State  Problem  in  One  Piinension 

In  the  previous  section  a general  development  of  the 
minimization  procedure  was  presented,  since  many  systems 
that  exhibit  non-adiabatic  effects  can  be  treated  as 
processes  occurring  on  two  electronic  potential  curves, 
a more  detailed  treatment  of  the  two-state  case  will  now  be 
given.  For  simplicity  it  will  be  restricted  to  one  dimen- 
sion, indicated  by  the  radial  variable  R. 

A general  form  for  the  real  unitary  transformation 
matrix  for  this  case  is  given  by 


(Kiel 

“(R|u) 


(2-39) 


with  v“(R^;a)  = 0,  y“(R^;o)  = j and  where  the  parameter  a is 
to  be  determined.  In  one  dimension,  the  momentum  operator 
is  given  by 


and  the  antisymmetric  matrix  of  the  adiabatic  momentum 
coupling  terms  has  the  form 


(2-41) 


T*(R|  . <*?|d*t/aB>  . (2-42) 

Using  Bgs . (2-39)  through  (2-42)  in  Eg.  (2-34)  leads  to 

j dR(T®(B)  ♦ t”(Rio)'1*  - 0 (2-43) 

where  the  "prime"  denotes  differentiation  with  respect  to  R. 
The  domain  of  integration  in  this  caee  will  be  the  interval 
(R^.R,!  (see  Fig.  (1()  where  t*(B)  is  significant.  In 
analogy  to  Eg.  (2-41)  one  has 


(2-31), 


Assuming  the  integrand  in  Eg.  (2-43) to  be  continuous  gives 
j dR[T®(R)  » y“ 


‘(R;n)'J3-r' 


A functional  form  of  y"  that  varies  smoothly  from  tero 
for  large  R to  r/2  for  small  R (see  the  discussion  in 
Section  2*1)  was  introduced  by  choosing 

Y“(R;a>  - (t/4)  11-tanhalR-Rjj)]  (2-47) 

where  R^^  is  the  point  where  the  diabetic  surfaces  cross. 


r"(Ria)'  = -(Xy'4)  aSech*o(R-Rjj)  (2-48) 


Sy"(P;o) '/8a 


(ir/2)  sech*a(R-8jj) 


• (a(R-Rjj)tanhQ(R-Rjj}-l/2J  . 


Using  Eqs.  (2-49)  and  (2-49)  1 


2-46)  leads  t 


dR  sech'a(R-Rjj)  (t®(R)  - (n/4)  sech^a  (R-R^,) ) 

Jt  |a(R-Rjj)tanha(R-8jj)  - 1/2]  = 0 . 


Equation  (2-50)  does  not  give  an  analytical  solution  for  the 
parameter  a ao  that  a numerical  procedure  must  be  implemented 
to  determine  the  value  of  a such  that  the  integral  is  less 
than  a small  positive  number. 


- 0 


Vl“(RI  - Vi*i(R)CQ8“f“lR,C)  + V,®, 
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v“a(R)  -Vj“i(Rl  = (V®i-Vj®  >ainY“lR,c)  •cosy“(R,c) 

where  (u,c|  equals  (m,a|  or  (d.R  ),  respectively. 

Equation (2-52)  emphasizes  the  point  tnade  in  the  last 
part  of  Section  2*2b  that  the  behavior  of  the  diabatlc 
potentials  depend  on  the  momentum  ooupling  i^Rl . 

Letting  R^  equal  R^  In  Fig.  (1),  one  sees  that  unless  the 
integral  in  Eq.  (2-52)  from  R^  to  R^  is  equal  to  it/2,  the 
diabatic  potentials  will  not  ooincide  with  adiabatic 
potentials  in  region  (a) . One  also  has  that  the  diabatic 
will  cross  whenever 

T‘*(R,Rg)  = l(2ntl)7i/4  , n = 0,1,2,...  (2-56) 

and  that  depending  on  t^r),  they  can  cross  more  than  once. 
These  undesirable  features  do  not  occur  in  this  procedure 
since  it  started  with  a proper  form  of  y®. 

2-5  A Model  Calculation 

Xn  the  following  calculations  t^IR)  was  given  the  form 
of  a Oaussian  and  written  explicitly  as 


Eq.  12-57) 


in  the  region  of  the  pseudocrossing.  Using 
Eq. (2-52) leads  to  the  result 


j/t/  (4b)  (erf  l.S(R-Rjj)  J-erf  IRjjSj^) ) ) 
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where  erf(x)  (Abramowitz  and  Stegun,  1972]  is  the  error 
function.  To  obtain  Che  parameter  in  Eq.  (2-50)  was 
integrated  numerically  using  Simpson's  rule (Abramowitz  and 
Stegun,  1972)  and  variable  step  sizes. 

Horse  potentials  (Horse,  1929),  defined  by 

V?i(R)  * D^(exp[-2a^(B-R°)  3-2  exp  (R-B°)  J ) 

1=1,2,  were  used  for  the  adiabatic  potential  curves  vfj  and 
V^j.  V^,  was  chosen  to  roughly  oorrespond  to  giving  a 
well  depth,  D^,  of  .176  a.u.  with  a minimum  located  at 
Ri  = 1.4  a.u.  and  a value  of  Si  of  .BOl  a.u.  vfs  was 
rather  arbitrarily  chosen  to  give  a well  defined  avoided 
croesing  at  a distance  of  2.2  a.u.  This  is  somewhat  like 
the  pseudocrossing  between  R*^-B2  and  h-Hj  (Tully  and 
Preston,  1971).  The  values  used  were  .139  a.u.  for  Dj, 

2.1  a.u.  for  Rf  and  1.9  a.u.  for  a,. 

The  value  of  b used  in  Eg.  (2-57)  for  t*(R)  was  50  a.u. 
and  was  set  equal  to  2.2  a.u.  Hith  this  choice  of  b,  the 
Gaussian  in  Eg.  (2-57)  has  a half  width  of  approximately 
0.2  a.u.  This  choice  of  b reflects  the  rather  sharp  avoided 
crossing  of  the  adiabatic  potentials.  Three  different 
values  of  were  used  in  this  study  all  with  b = 50  a.u. 

The  value  of  = 10  a.u.  roughly  corresponds  to  the  re- 
sults shown  in  Fig.  (2)  in  the  paper  by  Tully  and  Preston 
(Tully  and  Preston,  1971).  The  values  of  t equal  to  7.5  a.u. 
and  5 a.u,  were  included  for  comparieons. 


For  T equal  to  10  a.u.  the  value  of  the  parameter  a 
aatiafying  Eq.  (2-SO)  was  found  to  be  10. B9  a.u.  The  re- 
sults of  Eqs.  (2-S31  through  {2-551  for  uam  and  d are  shown 
in  Fig.  121.  One  notices  that  the  diagonal  curves,  and 
Viz  cross  twice  and  they  don*t  coincide  with  the  adiabatic 
curves  for  small  values  of  h.  Also  as  seen  in  Fig.  (3)  the 
off  diagonal  matrix  elements  Vi^z  diverge.  This  procedure 
clearly  avoids  these  disparities. 

The  iDOiQentum  couplings  t™  given  by  Eq.  (2-451  and  T® 
from  Eq.  <2-571  are  compared  in  Fig . (4|.  One  sees  that  this 
procedure  reduces  the  couplings  by  about  a factor  of  three. 
Also  shown  in  the  figure  is  the  absolute  value  of  the 
derivative  with  respect  to  H of  y“,  which  would  be  equal  to 
T in  the  strictly  diabatic  case.  From  Eq. (2-35 1 , the  range 
of  energy  where  these  couplings  can  be  neglected  is  given  by 

Ergi<<(M/21  |v“/P,“,  I . (2-601 

From  Fig.  (31  and  (41 , estimates  of  v"  and  Pi"  ace  given  by 
.002  a.u.  and  3 a.u.  using  that  the  reduced  mass  M of  the 
H -Hz  system  is  1836  a.u. , and  substituting  the  values  in  Eq. 
(2-601  justify  neglecting  the  momentum  coupling  if  is 
much  less  than  about  .0004  a.u.  (or  . 25  kcal/molel . Then  for 
very  low  collision  energies  our  procedure  would  be  useful  whereas 
for  higher  energies  an  adiabatic  representation  may  be  more 
convenient.  This  requirement  of  very  low  energies  is  due 
partially  to  the  small  reduced  mass  for  h'^-Hz  , and  should  not  be 
as  severe  for  other  systems  with  larger  reduced  masses. 


Fig.  {3J  A comparison  of  the  off  diagonal  potential  matrix 
elements  for  the  two  formulations  of  diabatic 
states  for  « 10  a.u.  The  dash-dot  and  dashed 
curves  correspond  to  and  VOj  in  Eq.  (55) 
respectively,  t*  given  by  Eq.  (57)  corresponds 
to  the  solid  line  and  is  included  for  comparison. 
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equal  to  7.5  a.u.  the  value  o£  the  parameter  a 
waa  found  to  be  9.172  a.u.  The  diagonal  elements  of  and 
Yni  compared  in  the  crossing  region  in  Fig,  (5). 

Although  the  striotly  diabatio  curves  do  not  have 
multiple  crossings  in  this  case,  they  still  don't  give  the 
correct  behavior  for  small  R.  As  shown  in  Fig.  (6),  the 
off  diagonal  elements  of  V again  diverge. 

The  momentum  couplings  i®  and  t”  are  compared  in 
Fig  (7).  one  notices  that  the  momentum  coupling  is  reduced 


e energy  r 

coupling  can  be  neglected  is  given  n 
(or  4.4  koal/mole>.  Thus  for  this  c 
has  considerably  increased. 

Figures  8) through  (101  show  the  results  for  i - 5 a.u. 
In  this  case  a is  6.840  a.u.  Pig,  (8)  again  demonstrates  the 
improper  behavior  for  small  R of  the  Vj^^'s.  One  also 
notices  that  as  is  decreased  the  crossing  from  this  pro- 
cedure occurs  over  a larger  region.  This  coincides  with 
the  larger  spread  in  v/°i  shown  in  Pig.  (9).  one  also 
notices  the  divergence  of  Vj'\  in  this  figure.  As  noted 
previously,  this  will  always  occur  if  the  integral  of  t®(r) 


The  coupling  terms  are  compared  in  Pig.  (loj  where  t 
sees  that  this  procedure  reduces  couplings  by  a factor  of 
around  ten . with  V ."j  and  P 
respectively,  the  energy  must 


(■'’■0)  A 
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Fig.  t9)  A comparison  of  the  off  diagonal  potential  matrix 
elements  for  the  two  fomulations  of  diabetic 
states  for  = 4 a.u.  The  dash-dot  and  dashed 
curves  correspond  to  VJj  and  in  Eq.  (S5) 
respectively,  t*  given  by  Eq.  (S7)  corresponds 
to  the  solid  line  and  is  included  for  comparison. 
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(or  19  )ical/mole).  Thus  again  the  range  in  energy  where  the 
momentum  couplinga  can  be  neglected  has  considerably 
increased. 

2-6  Discussion 

This  chapter  has  proposed  an  alternative  procedure  to 
construct  electronic  states.  This  procedure  is  based  on 
choosing  a transfortnation  from  the  adiabatic  representation 
into  another  that  will  ensure  the  proper  behavior  of  the 
potentials  away  from  the  crossing  region.  The  new  trans- 
formation depends  on  a set  of  parameters  which  are  obtained 
through  Eg.  (2-30).  It  was  shown  in  Section  2-4  that  this 
procedure  led  to  a rather  simple  treatment  of  the  two  state 
problem  in  one  dimension  while  the  results  demonstrated 
that  it  gives  well  behaved  potentials  and  reduced  momentum 
couplings. 

The  transformation  only  depends  on  the  momentum 
couplings  through  the  set  of  parameters  (a.).  The  problem 
is  thus  broken  up  into  two  parts:  (1)  the  determination  of 

the  parameters  (which  gives  also  the  momentum  couplings), 
snd  (2)  the  use  of  the  transformation  matrix  to  obtain  the 
diabatio  representation.  In  contrast,  in  the  strictly 
diabatic  transformation  a numerical  procedure  must  be 

X.  The  transformations  in  this  treatment  are  alwaya  analy- 
tical so  chat  numerical  procedures  need  only  be  used  once 
to  determine  the  parameters. 


As  was  emphasized  earlier,  whether  the  momentum 
oouplings  can  be  neglected  or  not  depends  on  the  system 
being  investigated  and  the  collision  energy.  An  appealing 
aspect  of  Che  present  procedure  is  chat  it  allows  us  to 
estimate  the  range  of  energy  where  it  can  he  used.  If  the 
collision  energy  is  not  in  this  range,  one  could  use 
either  the  adiabatic  representation  with  its  diagonal 
electronic  Hamiltonian  matrix  or  the  diabatic  representa- 
tion, Eq. (2-15),  with  the  possibility  of  unrealistic 


potentials. 

Several  related  problems  can  be  studied  along  the 
present  lines.  In  particular,  a numerical  treatment  of 
the  similar  problem  for  three  electronic  states  in  one 
dimension  would  be  useful.  Extensions  of  the  formalism 
would  also  be  helpful  in  describing  reactive  scattering. 
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and  the  notation  wae  previously  defined  in  Section  2-2. 
Since  there  is  only  one  electronic  state,  i.e.  no  coup- 
lings. the  adiabatic  representation  will  be  used.  Using 
the  fiorn-Oppenheimer  approxlmtion,  the  solutions  to  Eq. 
(3-1)  are  written  as 


I»(R)>  = |♦(H)>V(E^.  (3-3) 

Replacing  this  in  Eg.  (3-1),  multiplying  from  the  left 
with  <4(R)|  and  integrating  over  electronic  variables 


(3-4) 


V(R)  = <*(R)  |H^j^|*(R)>,  (3-5) 

i.e.  the  adiabatic  potential  energy  surface.  It  has  been 
assumed  that  the  electronic  functions  are  real  and  normalized. 
Equation  (3-4)  is  the  usual  time  independent  Schrodinger 
equation  for  one  electronic  state  systems  which  is  the 
starting  point  for  many  treatments  of  this  problem. 

Using  the  nuclear  coordinate  system  described  in 


Section  2-2,  Eq.  (3-4)  becomes 


SI 

- +V{R)(l(R)  = ElMR).  (3-61 

A form  for  the  solution  of  Eq.  (3-6)  is  appsrent  if  one 
considers  the  case  of  V(R)  equal  zero.  In  this  case  the 
solutions  are  just  traveling  waves.  A form  similar  to 
this  will  be  chosen  for  the  case  of  a nontero  potential. 
This  is  referred  to  as  the  Eikonal  Approximation  and  the 
solution  is  written  as 

^(R) 

*>(R)  - x<?)e  ' • (3-7) 

where  S(R)  Is  the  eikonal,  assumed  along  with  x(R)  to  be 
teal,  nstng  this  in  Eq.  (3-6)  leads  to,  after  some  sim- 
plification, 

■^(VS)-('7x)  - f5(v2s)*  + (V-E)x  = 0 

(3-8) 

where  the  R dependence  has  been  suppressed.  The  essential 
need  at  this  point  is  to  find  an  auxiliary  equation  in- 
volving the  eikonal  that  will  lead,  on  a transformation 
into  time,  to  a straightforward  classical  interpretation 
of  it.  An  expression  that  will  be  seen  to  fulfill  this 
is  obtained  by  requiring  that  the  gradient  of  the  ei)tonal 
satisfy  a Bamilton-Jacobi  like  equation,  i.e.  define 


(3-9) 


where  w is  to  be  determined.  Using  this  in  Eq.  (3-S) 


^ l’S)-l7x)  - ^ (V^Sjx  -f  IV  - M)x  - 0 .(3-10) 

An  expression  for  W is  obtained  from  the  real  part  of  Eq. 
(3-10)  and  one  finds  that 


(3-11) 


The  imaginary  part  of  Eq.  (3-10)  would  lead  to  an  expres- 
sion for  the  flux  but  since  it  is  not  necessary  for  this 
discussion,  the  treatment  of  it  will  be  deferred  to  a 


The  next  step  involves  implementing  the  short  wave- 
length approximation.  This  approximation  essentially 
assumes  that  the  nuclear  wave  functions,  Xi  vary  slowly 
so  that  terms  involving  oan  be  neglected.  It  can  be 
shown  (this  will  be  treated  in  detail  in  Section  3-S) 
that  if  this  approximation  is  used  in  Eq.  (3-10)  and  a 
transfomiation  into  time  is  made,  x will  satisfy  a first 
order  differential  equation  in  time.  Even  Chough  the 
case  being  considered  here  is  trivial  in  that  no  elec- 
tronic transitions  take  place,  similar  results,  i.e. 
first  order  differential  equations,  will  be  obtained  in 
t.he  general  case.  Further,  as  will  be  shown  in  Section 
3-6,  for  X expressed  in  several  convenient  forms,  these 


differential,  equations  will  result  in  x beinq  determined 
by  Hamilton’s  equations  of  motion. 

Using  the  short  wavelength  approximation  in  Eq.  (3-11) 
and  replacing  Che  result  in  Eq.  (3-9)  leads  to 


(3-121 


On  transforming  into  time,  the  gradients  of  the  eiXonal 
become  the  nuclear  momenta,  P(t),  and  the  nuclear  positions 
become  functions  of  time,  R(t),  so  that  Eq.  (3-12)  becomes 

♦ V(R(tl)  = E (3-13) 


Defining  the  Hamiltonian, 

H(R(t),  P(t))  = + V(H(t)), 


choosing  Che  initial  positions  R(t.  ) and  m 


(3-15) 


and  requiring  t 


1-16) 


i-e.  that  the  nuclear  poaitiona  and  cnomenta  aatiafyin^ 
Hamilton'a  equations  of  notion  result  in  conservation 
of  total  energy  (this  can  be  easily  seen  by  taking  the 
total  tine  derivative  of  the  Hamiltonian  and  using  Eqs. 


The  preceding  discussion  has  demonstrated  that,  for 
the  one-surface  case,  this  formalism  leads  tc  the  well 
known  classical  trajectory  method.  That  is,  the  nuclei  a 
treated  as  classical  particles  that  evolve  on  a quantum 
mechanical  potential  energy  surface.  This  should  aid  in 
seeing  through  some  of  the  complications  that  arise  in  On 
more  general  cases  to  be  considered  next. 

3-2  General  Time  Independent  Formalism 

In  the  preceding  section,  development  of  the  one- 
surface  case  was  presented  and  it  was  shown  that  the 
formalism  led  to  the  method  of  classical  trajectories. 

The  rest  of  this  chapter  w 
system  of  n electron  states, 
developed  without  reference  t 


1 develop  the  formalism  for  a 
The  formalism  will  be 
a particular  coordinate 


B generality  of  t 


system  in  order  to  emphasise 
dure  and  avoid  cumbersome 

In  order  to  avoid  assumptions 
of  the  nuclear  positions  and  momenta,  this  procedure  < 
nates  in  a time  independent  formalism.  Specif icially, 
seeks  the  solutions  of  th 
equation  which  satisfy 


time  dependence 


time  independent  Schrodinger 


[3-18) 
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(H-El |f>  - 0 

where  | > refers  to  electronic  coocdinates,  < | > indicates 
integration  over  electronic  coordinates. 


and  the  terms  in  Eq.  [3-19)  were  defined  in  Section  2-2. 

For  brevity,  the  dependence  on  nuclear  coordinates  has  been 
suppressed.  As  was  discussed  in  chapters  one  end  two,  one 
expands  the  solution  of  Eq.  (J-18)  in  an  electronic  basis. 
In  order  to  further  clarify  the  problem  of  non-Hermiticity 
mentioned  in  Sections  1-2  and  2-2c,  a specific  electronic 
representation  will  not  be  used  at  this  point.  The  basis 
will,  however,  be  assumed  to  be  real,  complete  and  ortho- 
normal. Expanding  in  this  arbitrary  electronic  basis  leads 


[»>  - !♦>  » (3-20) 

where  !♦>  is  a Ixn  row  matrix  and  i|r  is  a nxl  column  matrix. 
Using  Eg.  (3-20)  in  Eg.  (3-18),  multiplying  from  the  left  by 
<d|  and  integrating  over  electronic  coordinates  gives 


(3-22) 
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(3-23) 


and  one  nust  reioeinber  that,  due  to  the  presence  of  the 
gradient,  the  matrix  v must  be  treated  as  a vector.  Obvi- 
cuely  Eqs.  (3-221  through  (3-24)  define  nxn  square  matrices. 

As  was  mentioned  previously,  Eq.  (3-21)  is  not  in 
general  Hermltlan,  The  terms  that  give  rise  to  this  come 
from  Eq.  (3-22).  Since 


t difficult  to  s 
t)ie  electronic  b 


, by  making  use  of  the  complete- 


(3-27) 


£*''S  • ■ ! 
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This  equation  is  Hermitian  and  provides  the  starting  point 
for  the  developnents  to  be  presented  in  the  next  sections. 
3-3  The  Eikonal  Approximation 

As  in  the  one  surface  case,  the  next  step  in  this 
treatment  involves  using  the  Sihonal  Approxiroation.  The 
form  of  the  solution  is  chosen  to  be  the  product  of  an 
amplitude  written  as  a nxl  column  matrix  of  complex  ele- 
ments and  a common  phase  that  is  proportional  to  the 
eikonal.  A common  phase  is  used  because  it  will  lead  to 
eoraraon  momenta  and  positions  for  all  electronic  channels. 
Explicitly,  the  form  of  the  nuclear  wave  function  is  chosen 

J(R)  - J(?)e”  " (3-31) 

where  S(R)  is  the  eikonal.  Using  Bq.  (3-31)  in  Eg.  {3-30} 


(3-32) 

♦ (7S)  -vx  ♦ j v*?x  + § v*vy  + Vx  - Ey  = 0 

As  in  Section  3-1,  the  gradient  of  the  eikonal  la  required 
to  satisfy  a Kamilton-Jacobi  like  equation  which  gives 


3-33) 


S9 

where  W is  to  be  determined.  Using  Eg.  (3-33)  in  Eg.  (3-32) 

- S 

(3-34) 

♦ & v-7x  + I v-vx  * J'J  ■ "J  ■ 0 • 

An  expression  for  if  can  be  obtained  by  multiplying  Eg.  (3-34) 
from  the  left  by  multiplying  the  adjoint  of  Eg.  (3-34) 
from  the  right  by  x and  adding  the  results,  Doing  this 

W = y'X^X)*^(2^^Vj  + ^ 7S- [(7x^)X  - X^OX) 

' + I'/v-vx  - (''x''')-vx) 

* ""J  * 2’®'jS^;:x1  • (3-3S) 


This  expression  for  H may  seem  to  be  a bit  complicated  and 
perhaps  a few  comments  on  its  general  properties  may  be 
helpful, 

considering  an  initial  channel  where  one  of  the 
amplitudes  eguals  one  and  all  others  are  zero,  then  if  an 
initial  time  is  chosen  so  that  the  interactions  between 
the  fragments  and  the  couplings  (v)  are  zero  and  the 
amplitudes  are  constant,  the  expression  for  w reduces  to  a 
single  potential  energy  surface.  This  is  certainly  what 


one  would  expect.  As  time  increases  and  the  fragments 
enter  the  interaction  zone,  the  amplitudes  are  free  to 
change.  In  this  case  there  are  three  contributions  to  S. 

The  first  contribution  comas  from  the  first  term  on  the 
right  hand  side  of  Eq.  (3-35).  This  can  be  viewed  as  an 
average  potential  with  the  weights  determined  from  the 
amplitudes.  The  second  contribution  comes  from  the  next 
two  terms  in  Eq.  (3-35)  and  arises  from  the  variations  of  x. 
If  the  elements  of  j are  slowly  varying,  this  second  contri- 
bution may  be  'small*  enough  to  neglect.  The  final 
contribution  comes  from  the  last  three  terms  which  depend  on 
the  electronic  couplings.  At  a final  time  when  the  ampli- 
tudes are  again  constant  and  the  electronic  couplinge  and 
interactions  between  the  fragments  are  zero,  <7  Is  again  given 
by  the  first  tern  on  the  tight  hand  side  of  Eq.  (3-35).  How- 
ever w for  the  final  state  is  now  an  average  potential  with 
the  weights  detemlned  by  the  final  values  of  the  amplitudes. 
The  implications  of  this  will  be  discussed  in  a later  section 
after  a transformation  to  time  has  been  made. 

An  expression  for  the  flux  is  obtained  by  multiplying 
Eq.  (3-34)  from  the  left  by  , multiplying  the  adjoint  of 
Eq.  (3-34)  from  the  right  by  X and  ai;btracting  the  results. 
Carrying  out  the  algebra  leads  to 

- 2^?’x!  * = o . (3-36) 


Defining 


i,  = 5E 


and  recalling  that  the  electronic  basis  aatiefles  Eq. 
(3-23)  one  obtains 


(3-41) 


From  Eq.  (3-40)  one  sees  that  there  are  three  terms  in  the 
expression  for  flux.  The  first  tern  J defined  by  Eq. 

(3-37)  is  the  most  transparent.  At  the  initial  time  when 
V IS  the  zero  matrix  and  is  constant.  Bq.  (3-40)  reduces 
to  the  gradient  of  equal  to  zero,  j is  just  the  produot 
* velocity  term  given  by  (7S)/»  times  the  number  density 
particles  given  by  ^ x-  ht  later  times  when  v is  non- 
zero and  X is  changing,  the  flux  has  additional  terms 


involving  electronic  couplings  and  variations  of  If  x 
is  slowly  varying,  it  may  be  possible  to  neglect  the  last 
term  in  Eg.  (3-40)  so  that  the  flux  will  be  given  by  Eg. 

As  has  been  suggested  throughout  this  section,  if  use 
is  made  of  the  slowly  varying  nature  of  x,  the  expressions 
for  W and  J are  considerably  simplified.  This  will  be  the 
topic  of  the  next  section. 

3-4  The  Short  Wavelength  Maroximatlon 

In  the  previous  section,  the  Eikonsl  approximation 
was  used  and  it  was  found  that  the  nuclear  wave  function, 

X,  satisfied  Eg.  [3-34).  This  section  will  make  use  of  the 
slowly  varying  nature  of  x in  order  to  introduce  a physi- 
cally motivated  approximation  to  Eg.  (3-34).  This  will  in 
turn  simplify  the  expressions  for  if  and  the  flux. 

The  approximation  referred  to  above  is  the  well  known 
short  wavelength  approximation  [Newton,  1966).  This  approx- 
imation is  based  on  the  classical  behavior  of  the  nuclear 
particles.  From  Eg.  [3-33)  one  has 

|9S|  = /Zm(E-B)  (3-42) 

so  that  a characteristic  wavelength  can  be  defined  es 


a full  classical  treatment  of  the  trajectories,  fi  would 


the  wavelength  would 


should.  The  short  wavelength  approximation  assumes  that 


which  is  normally  valid  except  around  turning  points  where 
E » W.  Obviously,  from  Eq.  (3-43),  as  the  total  energy  is 
Increased  this  approximation  becomes  better  so  that  it  is 
in  actuality  a high  energy  approximation. 

In  a similar  fashion  one  can  define  "characteristic 
wavelengths’  for  the  amplitudes  as 


a. 


for  an  arbitrary  electronic  channel  j.  If  the  amplitudes 
are  assumed  to  be  slowly  varying  functions,  their  "wave- 
lengths" should  be  large  compared  to  I so  that  with  Eq. 
13-44)  one  would  have  that 


(3-46) 


(3-47) 


*j,2  * 

using  these  conditions  in  Eq.  (3-34)  leads  to 

- ^(7S)  • (7x>  - t (7S)  -n  ♦ t yvx  + Vx  - HX  = 0 . 

(3-46) 

Carrying  through  the  same  procedure  that  was  used  to  obtain 
Eq.  (3-3S)  results  in 

W = + ^(7S).((7x'^)X  - /vx) 

+ 2(7S)-j'^vx  + mx^v-Vj()  . (3-49) 

Assuming  that  the  second  term  on  the  right  hand  side  of 
Eq.  (3-49)  is  negligible  gives 

w - (x^|)‘V‘j  * * I r?’ 2 • 

The  general  properties  of  H discussed  in  the  previous 
section  still  hold  here.  Suffice  it  to  say  that  the 
approximations  used  to  obtain  Eq.  (3-SO)  have  considerably 
simplified  the  determination  of  w in  that  only  a )cnowledge 
of  the  amplitudes  and  not  their  variations  is  needed. 

The  flux  is  determined  by  using  the  same  procedure 
that  led  to  Bq.  (3-36).  Carrying  through  the  algebra 


5-2,  = 0 (3-52) 

where  is  given  by  Eq.  (3-37).  Equation  (3-S2)  is  an  expres- 
sion for  the  conservation  of  current.  In  the  case  of  one 
electronic  surface,  Bq.  (3-37)  reduces  to 


2,  = X*  ^ ' '3-53) 

which  is  a well  known  result  (Messiah,  1966). 

Equations  (3-48) , (3-50)  and  (3-52)  are  the  main  results 
of  this  section.  One  notices  that  Eq.  (3-48)  does  not  con- 

section  that,  because  of  t)iiB,  the  amplitudes  will  satisfy 
first  order  differential  equations  in  time.  This  is  a 
pleasing  aspect  since  Hamilton's  equations  of  motion  are 
also  first  order  differential  equations  in  time. 

Another  point  worth  mentioning  has  to  do  with  the 
form  of  if.  At  a later  stage  in  this  development,  if  will 
be  used  to  construct  a "classical"  Hamiltonian  and  one 
sees  from  Eq.  (3-50)  that  W depends  on  the  amplitudes. 


Then  the  ■classical"  Hamiltonian  will  not  only  depend  on 
the  traditional  variables  hut  also  on  the  nuolear  amplitudes 
as  well.  This  suggests,  already  at  this  point,  that 
all  the  nuclear  dynamics  could  perhaps  be  determined 
through  Hamilton's  equations  of  motion. 

3-5  Time  Peocndent  Equations 

The  developments  of  the  preceding  sections  of  this 
chapter  were  done  with  a time  independent  formalism.  This 
section  shall  be  concerned  with  transforming  these  previous 
results  into  a time  dependent  picture.  The  use  of  the 
word  time  is  somewhat  arbitrary  but  it  will  be  seen  later 
that  the  parameter,  which  the  nuclear  variables  are  a 
function  of,  plays  the  role  of  a *claaslcal"  time. 

In  making  a transformation  into  time  one  has 

a -■  R(t)  , (3-54) 

VS  - 7S(t)  = P(t)  (3-55) 


j(B(t))  A(t)  . (3-55) 

Matters  are  somewhat  simplified  if  a coordinate  system  is 
chosen  such  that  one  of  the  orthogonal  unit  vectors,  say 
s,  is  in  the  direction  of  P (see  Fig.  11).  Then  one  has 


p = IPI 


(3-58) 


7's  = n l^dnP  ) 


(3-59) 


PurthiernLore, 


^ u)»  u V . as  at  _ - 

3t3?i”dt3s3t  3% 


(3-60) 


so  that  Eq.  (3-48)  becomes 


(3-61) 


By  the  now  standard  technique,  one  can  show  that 


(3-62) 


A somewhat  simpler  differential  equation  that  Eq.  (3-61) 
results  if  one  defines 


A(t)  -C(t)exp^(  I dt'(W(t‘)  -Sj^.lnp'')) 
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where  t,  is  the  Initial  time,  so 


T * '^  * !•?  * I “ 2 • 

It  is  also  relatively  easy  to  show  that 

W = (C*C)*V{V  + P.v  + I vv)C  . (3-65) 

By  multiplying  Eq.  (3-64)  from  the  left  by  C*  and  the 
adjoint  of  Eg.  (3-64)  from  the  right  by  C and  adding  the 
results,  one  finds  that 

^(C^C)  = 0 (3-66) 

i.e.  the  probability  is  conserved.  Equation  (3-63)  can  also 


A(t) 


(3-67) 


which  shows  that  the  amplitudes  have  the  correct  asymptotic 
form.  The  amplitudes  C are  more  convenient  to  wor)c  with 
because  they  avoid  the  singularities  that  occur  in  A at 
the  turning  points  where  P is  zero. 

The  main  result  thus  far  is  given  by  Sq.  (3-64).  If 
the  trajectories  and  momenta  are  Itnown  functions  of  time. 


1 )3-3) . reqi 
Bq.  (3-33)  1 


- ii  • 


(3-71) 


s-f  . 

s=-#  , 


I (3-72)  and  (3-73) 


(3-75) 
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chcsen  to  be  conjugate  variables,  then  with  trajectories 
and  momenta  satisfying  equations  13-72)  and  (3-73)  re- 
spectively, energy  is  conserved. 

It  is  not  difficult  to  show  that  for  the  amplitudes 
expressed  in  terms  of  their  real  and  imaginary  parts 
Egs.  0-74)  and  (3-75)  are  equivalent  to  Eg.  (3-54).  The 
advantage  offered  by  Eq.  (3-«4)  is  that  the  fora  of  the 
amplitudes  is  not  specified.  It  will  be  shown  later  that 
the  flexibility  in  choosing  the  form  of  the  amplitudes  will 
lead  to  other  sets  of  canonical  variables  that  have  a 
straightforward  physical  interpretation. 

The  next  section  shall  be  concerned  with  finding 
solutions  of  £q.  (3-64).  Several  different  forms  for  Che 
solution  will  be  developed  and  it  will  be  seen  that  they 
satisfy  Eqs.  (3-74)  and  (3-75).  Having  more  than  one 
form  is  certainly  not  surprising.  It  is  in  a way  equivalent 
to  the  problem  in  classical  mechanics  of  having  more  than 
one  choice  of  coordinates  some  of  which  lead  to  a simpler 
treatment  of  the  problem. 

3-6  Solutions  of  the  Equations 

This  section  will  develop  solutions  to  Eq.  (3-64)  for 
two  forms  of  the  amplitudes.  It  will  be  shown  how  these 
forms  can  be  used  to  construct  canonical  variables  that 
satisfy  Eqs.  (3-741  and  (3-75).  It  will  also  be  emphasised 
that  care  must  be  used  in  selecting  a set  of  canonical 
variables  in  that  some  of  them  can  lead  to  inherent  diffi- 
culties that  arise  from  the  form  of  the  Hamiltonian. 


in  general  a mcl  column  matrix  of 


The  matrix  C is 
complex  numbers.  A general  way  of  writing  such  a matrix  is 
to  let 


C{t) 


c(t) 


(3-76) 


where  y Is  a real  square  diagonal  nxn  matrix  and  c is  a 
real  nxl  column  matrix.  Using  Eq.  (3-76)  in  Eq.  (3-64)  and 
rearranging  terns  give 


s j " component  satisfies 


3 ft  ■ fi  ® 


An  equation  for  o^  can  be  found  by  multiplying  Eq.  (3-7B) 
from  the  left  by  c^,  multiplying  the  complex  conjugate  of 
Eq.  (3-78)  from  the  left  by  and  adding  the  results. 
Doing  this  results  in 


(3-79) 


• S ®3®*  ^ Ji  "jk®  *■  ''  ®k 

”3k  ' ’'jk  + £*''jk  + I I ’'jj*  ''jit  • [3-eo) 

With  the  form  of  C given  by  Eq.  (3-76),  one  has  that 

w - (c^c)*^  o^ e^*  S ® ' (3-81) 

where  the  elements  of  U were  defined  in  Eq.  (3-80).  It  is 
straightforward  to  show  that 

iS.=  (T.  o »«’^3  r n *5’'^ 

TyJ"  f=<=jl  I O.^e  o^. 


Comparing  Eqs.  (3-79)  and  (3-82)  gives 
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or  by  defining 


end  recalling  from  Eq.  0-66)  that 


1S_ 


0-86) 


One  would  expect  at  this  point  that  and  are  canonical 
variables  and  that  they  should  satisfy  equations  similar  to 
t3-74)  and  0-75)  respectively.  The  other  relation  can  be 
obtained  by  sutractlng  the  results  that  were  used  to  obtain 
(3-79) , This  results  in  a differential  equation  for 
Yj  and  one  finds  that 


O-i 


(3-84) 


W as  a function  of  and  Yj  becomes 


* i T n p 

2 i ® 
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Comparing  equations  (3-87)  and  (3-90)  gives  the  desired 
result  in  that  one  has 


(3-91) 


Then  as  was  expected,  Yj  and  form  a set  of  canonical 
variables.  Expressing  the  total  time  derivative  of  the 
Hamiltonian  in  this  set  of  variables  and  using  Eqs.  (3-86) 
and  (3-91)  lead  to  energy  conservation. 

An  Interesting  aspect  of  these  variables  is  seen  by 
considering  Eq.  (3-69).  Since  the  terms  on  the  right  hand 
aide  contain  the  quantity  Pj*®, there  is  an  essential  singu- 
larity in  this  equation  if  P^  is  equal  to  zero.  This  can 
he  very  troublesome  since  normally  one  starts  in  a pure 
electronic  state  where  all  probabilities,  are  zero 
except  one  of  them  which  is  set  equal  to  one.  Trying  to 
integrate  this  set  of  coupled  equations  that  contain 
singularities  at  the  boundary  conditions  is  clearly  not 
meaningful. 

Even  though  this  set  of  variables  has  . the  difficulty 
that  was  mentioned  above  they  ace  easily  underetood  in  a 
physical  sense.  Indeed,  it  is  gratifying  to  see  that  the 
quantum  mechanical  probability,  p^ , and  phase,  Yj.  turn 
out  to  be  canonical  variables  In  a "classical'  description. 
Although  other  sets  of  canonical  variables  which  avoid 
these  difficulties  are  used  to  do  the  calculations,  it  is 


convenient  to  construct  the  probabilities  and  phases  as  a 
guide  to  a straightforward  physical  interpretation  of  the 
collision  process. 

Another  form  for  the  amplitudes  is  obtained  by  letting 


(3-92) 


Using  this  form  of  the  amplitude  in  Eg.  (3-64)  gives 
oh 

X at  * 3t  * “ 2 (3-93) 

where  the  elements  of  U axe  defined  by  Eg.  (3-80)  and  X 
and  Y are  real  nxl  column  matrices.  It  is  straightforward 
to  show  that  the  component  satisfies 

nh 

xar  * at^*  “pjh\  * IBl  Vk  “ “ • 

A differential  eguation  for  X.  is  obtained  by  subtracting 
Eg.  (3-94)  from  its  complex  conjugate  which  gives 

Jt*  ■ ra:  1“!  I”j,  - v»i.  • i«  . Oj’jlS,,  . I3-.5I 

In  a similar  fashion,  the  differential  eguation  for  is 
found  by  adding  Eg.  (3-94)  to  its  complex  conjugate  and 
one  finds  that 


..  (3-96) 


dY.  , . 

3i^  ■ - iw  >■!  • "ji.i’S,  • “ ! »»  - • 

Eijuationa  (3-95)  and  (3-96)  provide  first  order  differential 
equations  for  determining  the  amplitudes  but  the  taa)c  still 
remains  to  determine  sets  of  canonical  variables.  In  order 
to  find  sets  of  canonical  variables,  it  is  necessary  to  find 
an  expression  for  W.  to  somewhat  simplify  matters,  Bq. 

(3-66)  will  be  used  and  it  will  be  assumed  that 


The  results  of  this  analysis  do  not  depend  on  this  assump- 
tion in  that  if  Bq.  (3-97)  was  not  satisfied,  then  new 
variables  weighted  by  the  inverse  of  the  square  toot  of 
Bq.  (3-97)  oould  be  formed  that  would  also  satisfy  Eqs. 
(3-95)  and  (3-96).  Using  Eqs.  (3-9?)  and  (3-92)  in  Eq. 
(3-65)  gives 


« ' I • '3-98) 


It  is  not  difficult 


that 


w-  • -i»t"  I i«ii.  - 


I ‘“jk  * “jk>\  + I - “ik’^k  '3-1001 


where  use  was  nade  of  the 


property  that 


(3-101) 


Equations  (3-99)  and  (3-100)  are  alnilar  to  Eqs.  (3-95)  and 
(3-96)  respectively  and  for  certain  choices  of  o and  6, 

Xj  and  Yj  will  play  the  role  of  canonical  variables. 

By  cooiparinq  Eq.  (3-95)  to  Bq.  (3-99)  and  Eq.  (3-96) 
to  Eq.  (3-100)  one  finds  that  and  Y^  will  form  a canonical 
set  of  variables  if 


Two  obvious  choices  of  a and  6 that  satisfy  Eq.  (3-102) 
are  given  by  letting 


(3-102) 


(3-103) 


(3-104) 


only  possibilities.  The 


where  a and  8 are  given  by  Eq.  (3-103)  will  be  referred  to 
as  the  symmetric  form  and  this  form  was  used  in  doing  the 
calculations. 

One  notices  that  Eqs.  (3-991  and  (3-100)  do  not 
contain  singularities  in  the  variables  and  Xj  so  that 
the  difficulties  that  occurred  when  using  the  probability 
and  phase  as  canonical  variables  are  avoided.  As  was 
mentioned  earlier,  the  probability  end  phase  do  have  the 
advantage  of  being  easy  to  physically  interpret  and  of 
being  formed  quite  easily  by  using  the  relations 

Pj  = (3-106) 


1 

Yj  - - fi  tan  (jjji)  . (3-107) 

The  solution  of  the  problem  will  then  consist  of  integrating 
the  coupled  equations  given  by  the  expressions  (3-72)  and 
(3-73)  and  the  corresponding  ones  for  the  amplitudes.  The 
amplitudes  will  be  chosen  to  have  the  form  of  Eq.  (3-92). 

The  probabilities  and  phases  are  found  by  using  Eqs.  (3-106) 
and  (3-107). 

3-7  Expressions  for  Observablea 

As  was  pointed  out  in  Section  3-3,  one  of  the 
curious  features  of  this  formalism  is  that  the  trajectories 
end  up  on  an  average  potential  energy  surface  with  the 
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wei9t>ts  determined  by  the  probabilities  . This  may  seem 
to  be  a serious  drawback  because  experimentally,  the  sys- 
tem begins  and  ends  up  in  a definite  electronic  state. 

This  feature  is  not  unique  to  this  formalism  in  that  this 
problem  also  occurs  in  other  semi-classical  treatments  of 
multi -surf ace  systems. 

Trajectories  that  do  go  from  one  surface  to  another 
can  be  found  by  considering  the  quantum  mechanical  ex- 
pressions for  the  transition  matrix  or  scattering  matrix 
which  can  be  expressed  in  terms  of  the  Initial  and  final 
states  of  the  system.  The  problem  then  consists  of  finding 
trajectories  that  go  smoothly  from  the  initial  to  the 
final  state.  Several  procedures  have  bean  developed  tor 
finding  these  trajectories  but  they  are  quite  complicated 
and  not  easy  to  implement,  h further  discussion  of  this 
problem  will  be  presented  in  chapter  six. 

The  approach  to  be  used  in  this  work  is  based  on 
accepting  average  potentials  for  the  trajectories  and  on 
using  the  's  to  construct  a total  electronic  transition 
probability.  The  problem  is  that  there  are  many  "classi- 
cal" trajectories  that  correspond  to  the  sane  quantum 
mechanical  initial  state.  One  would  then  expect  that  some 
sort  of  an  average  of  these  "classical"  probabilities,  Py 
would  correspond  to  the  quantum  mechanical  probability. 

This  becomes  more  clear  if  one  considers  the  case  of 
acoLlinear  collision  between  an  atom  and  a diatomic 
molecule.  The  "classical"  initial  state  would  consist  of 


specifying  a total  energy,  E,  which  is  the  sum  o£  the 
vibrational  energy,  of  the  diatomic  and  the  relative 
hinetic  energy,  the  distance  between  the  atom  and  the 
center  of  mass  of  the  diatomic,  the  vibrational  coordinate 
of  the  diatomic,  the  direction  of  the  diatomic  momentum 
and  the  initial  amplitudes,  fhe  guantum  mechanical  Ini- 
tial state,  however,  would  represent  the  diatonic  by  a 
wave  function,  i(p^,  corresponding  to  the  n^  vibrational 
The  total  electronic  transition  probability  can 

^j-*j'E.En>  ‘ J (3-1081 

where  X is  the  vibrational  coordinate.  Equation  (3-108) 
provides  a prescription  for  determining  the  total  elec- 
tronic transition  probability.  Final  probabilities  for 
different  values  of  X are  calculated  and  the  total 
probability  is  obtained  by  averaging  the  final  probabili- 
ties with  the  weights  given  by  Eg.  (3-108). 

The  procedure  for  determining  total  probabilities  in 
the  general  case  is  essentially  the  same.  For  variables 
that  are  treated  quantum  mechanically  by  probability  dis- 
tributions, one  determines  the  final  probabilities  for  a 
number  of  values  of  the  variable  and  averages  them  with  the 
weights  dependent  on  the  probability  distribution. 


CHAPTER  4 

ELECTRONIC  PROBLEM  FOB  Hj"* 
HE  ADIABATIC  REPRESENTATION 


molecule,  the  surface  for  a triatomic  system  would  require 
in  the  neighborhood  of  a thousand  nuclear  configurations. 


Even  if  the  electronic  calculation  were  to  only  take  ten 
seconds  per  configuration,  the  surface  would  still  require 
nearly  three  hours  of  calculation  time.  These  limitations 
are  inherent  in  the  ab  initio  approach  and  from  the  second 
limitation  alone  it  is  unlikely  that  this  approach  would  be 
feasible  for  obtaining  the  surfaces. 

At  the  other  extreme  of  the  spectrum  lie  empirical 
approaches  to  the  problem.  These  approaches  are  baaed  on 
using  arbitrary  functions  to  construct  a potential  energy 
surface.  The  function's  parameters  are  adjusted  until  the 
surface  passes  through  either  experimentally  known  or 
theoretically  calculated  points.  Even  though  these  ap- 
proaches bypass  the  need  to  do  a quantum  mechanical 
calculation,  quite  a bit  of  guesswork  is  involved.  If 
these  approaches  were  used  on  systems  that  didn't  contain 
a wealth  of  experimental  or  theoretical  information  there 
would  exist  a significant  risk  of  obtaining  unreliable 

The  final  clasa  of  approaches  to  be  considered  here 
are  the  semi-empirical  ones.  These  approaches  have  the 
advantage  of  starting  with  the  time  independent  Schrodlnger 
equation.  This  can  not  only  lead  to  a better  understanding 
of  the  approximations  used  but  also  offer  some  insight  into 
its  limitations.  The  number  of  seml-erapirioal  approaches 
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is  staggering  and  in  a way  reflects  the  great  diversity 
of  physical  syatens  being  studied. 

As  is  well  known,  most  of  the  effort  in  an  electronic 
calculation  is  spent  in  evaluating  the  two  canter  electron 
integrals.  The  contrion  trait  of  semi -empirical  approaches 
is  that  they  either  simplify  or  reduce  the  number  of  those 
integrals.  A well  known  approaoh  that  simplifies  the 
exchange  integral  is  the  Xa  method  (Slater,  1971).  There 
are  a number  of  other  methods  such  as  CNDO  that  simplify 
the  problem  by  reducing  the  number  of  electron  integrals. 

One  of  the  possible  difficulties  that  arise  in  using 
most  semi-empirical  methods  is  that  they  usually  employ 
only  one  electronic  configuration,  i.e.  a single  antisym- 
“etrized  product  of  molecular  orbitals.  This  may  give  a 
good  description  of  the  surface  tor  some  regions  of  the 
internuclear  coordinates  hut  usually  does  not  adequately 
describe  the  entire  surface,  which  is  necessary  for  the 
scattering  calculation.  Although  there  are  semi-empirical 
methods  that  enploy  configuration  interactions,  the  number 
of  internuclear  configurations  needed  to  determine  the 
surface  would  more  than  lUcely  make  using  these  procedures 

The  method  that  will  be  employed  here  is  called 
Dlatomica  in  Molecules.  It  was  introduced  by  Ellison 
(Ellison,  1963)  and  later  generalized  to  include  direc- 
tional bonding  (Kuntz  and  Roach,  1972j  Tully,  1973).  it 
has  fairly  recently  been  used  to  construct  a number  of 
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potential  surfaces  for  triatomic  systems.  Two  of  the  many 
examples  are  the  surfaces  for  LiHj  and  FBj  (Tully,  1973a|. 

This  method  utilises  the  fact  that  the  electronic 
Hamiltonian  can  be  written  in  terms  of  atomic  and  diatomic 
Hamiltonians.  This  leads,  as  will  be  seen  in  the  next 
section,  to  a solution  which  is  expressed  in  terms  of  an 
overlap  matrix  and  atomic  and  diatomic  energies.  The 
ground  and  excited  energy  levels  of  the  diatomics  are 
obtained  through  either  experiment  or  theory,  whether  or 
not  the  electronic  problem  can  be  solved  with  this  approach 
depends  on  the  availability  of  the  diatomic  energy  levels. 

Another  advantage  of  this  method  comes  from  the  fora 
of  the  electronic  basis  functions.  Since  the  basis  func- 
tions axe  expressed  as  an  antisymmetrised  product  of  atomic 
orbitals,  a little  manipulation  of  the  antisymmetriser 
leads  to  a valence  bond  description  of  the  electronic  prcb- 
lem.  Thus,  this  method  has  built  into  it  the  correct 
electronic  description  of  the  reactants  and  products. 

In  the  following  section  the  formalism  for  Diatomics 
in  Molecules  will  be  briefly  developed.  This  formalism 
will  be  applied  to  the  Hg*  system  and  the  eigenvalues  and 
eigenfunctions  will  be  the  topic  of  Sections  4-3  and  4-4 
respectively.  Section  4-5  will  focus  on  the  non-adiabatic 
couplings  which  are  a convenient  by-product  of  this  method. 
4-2  Method  of  Diatomics  in  Molecules 

The  goal  of  Diatomics  in  Molecules  (DIM!  is  to  make 
use  of  the  energies  of  the  diatomics  that  comprise  the 
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polyatomic  to  construct  the  polyatomic  energies.  If  the 
diatomic  energies  axe  all  known  experimentally,  this 
procedure  eliminates  the  necessity  of  doing  an  electronic 
calculation  except  perhaps  for  overlaps.  In  most  applica- 
tions the  Intesnuclear  dependence  of  the  overlaps  is 
neglected  sc  that  the  electronic  calculation  is  avoided 
altogether  and  the  electronic  problem  is  solved  with  only 
a knowledge  of  the  diatomic  energies. 

Since  there  are  already  a number  of  detailed  treatments 
of  this  subject  available  in  the  literature  (Ellison,  1963j 
Tully,  1973b) , only  a brief  outline  of  the  method  will  be 
Presented  here.  This  will  hopefully  )\ave  the  advantage  of 
displaying  the  main  conclusions  without  resorting  to  the 
rather  cumbersome  formalism  that  is  inherent  in  more 
detailed  treatments.  The  development  presented  hare  will 
somewhat  follow  the  one  given  by  Tully  (Tully,  1980). 

The  basis  functions  in  this  procedure  are  chosen  to  be 
antisymmetrised  products  of  atomic  functions.  For  a system 
ooraposed  of  n electona  and  N nuclei  where  at  least  three 
nuclei  are  aesumed  to  be  present,  the  basis  functions  can 


«„(l,---,n)  = Aj,»^(l,...,n)  (4-1) 

where  is  the  n electron  antisymmetriser  and  the  « are 
defined  as  (Moffitt,  1951) 
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■■■''  m 1”  - "n  + I---''")  («-2) 

where  the  etomic  functions  n are  assumed  to  be  antisymmetric. 
The  notation  has  been  somewhat  changed  here  in  order  to 
emphasize  the  importance  of  assigning  electrons  to 
particular  nuclei.  The  bracket  will  still  be  used  to 
indicate  integration  over  the  relevant  electronic  coordi- 
nates. Expanding  the  total  wave  function  in  terms  of  the 
basis  of  Eq.  (4-1)  gives 

y u,---,n)  - : (<-3) 

where  are  the  expansion  coefficients,  using  Eg.  (4-3) 
in  the  time  independent  Schrbdinger  equation  leads  to 

Hr  - SrE  (4-4) 


e adiabatic  energies  a 


Equation  (4-5)  wag  obtained  by  using  the  property  that  the 
antigymmetrizer  commutes  with  the  Hamiltonian. 

The  next  step  in  this  procedure  is  to  partition  the 
Hamiltonian  into  components  that  are  atomic  and  diatomic 
in  nature.  This  can  be  aocomplished  by  writing  (Ellison, 
1563) 


which  is  comprised  of  the  nuclei  K and  L and  the  electrons 
assigned  to  those  nuclei.  A similar  definition  holds  for 


operator  H has  the  following  matrix  representation 


(N-2) 


(4-7) 


where  H*^’  is  the  Hamiltonian  for  the  isolated  diatonic 


H***.  In  the  electronic  basis  defined  by  Eq.  (4-1),  the 


I f h'"’’ 


(4-8) 


(4-10) 


Since  the  total  antisymnetrizer  does  not  coiiuBute  with  the 
iragment  Hamiltonians,  and  H***,  Eg.  (4-6)  is  not 

term  by  term  Hermitian.  The  sura  will  however  remain 
Kermitian  if  a complete  set  of  electronic  states  are  used. 

Equation  (4-B)  is  of  utmost  importance  in  this  theory. 
One  sees  from  Eqs.  (4-9)  and  (4-10)  that  the  total  Hamilto- 
nian is  now  expressed  in  terms  of  fragment  Hamiltonians 
whose  operators  correspond  to  isolated  atomic  and  diatomic 
systems.  The  solution  to  the  problem  will  consist  of  ma)ting 
transformations  of  the  electronic  basis  given  by  Eg.  (4-1) 
to  bases  that  diagonalize  the  fragment  Hamiltonians,  i.e. 
whose  eigenvalues  correspond  to  the  experimental  or 
theoretical  energies. 

In  order  to  show  how  these  transformations  are  done,  a 
brief  development  for  the  operator  is  presented  here. 

For  clarity,  the  basis  functions  are  written  in  the  fora 


where  is  the  produot  of  all  atomic  functions  not 

centered  on  atom  K and  is  the  atomic  function  centered 

on  atom  K.  of  course,  the  electrons  are  still  assigned 
according  to  Eg.  (4-1).  Assuming  that  the  electronic 


where  the  4 are  the  expansion  coefficients.  Substituting 
Eq.  (4-121  Into  Sq.  14-9)  gives 


where  the  elements  of  S were  defined  in  Eq.  (4-5).  The 
important  thing  to  reoognire  in  Bq.  (4-13)  is  that  the 
matrix  4^*^^  depends  only  on  the  atom  K.  Then  it  should  be 
possible  to  relate  4**’  to  the  eigenvalue  problem  of 
for  isolated  atom  K lay  means  of  a transformation.  Since 
the  atomic  basis  is  assumed  to  be  complete  and 
operates  only  on  atom  X one  also  has  that 


(4-14) 


Defining 


h(K) 


(4-lS) 


one  obtains  the  matrix  equation 
h'*’  ■ 


»-17) 


Considering 
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the  eigenvalue  problem  for  isolated  atom  K, 


where  and  are  the  matrices  of  eigenvectors  and 
eigenvalues  respectively  and  x**'*  is  a row  matrix.  The  two 
basis  can  be  related  by  means  of  a transformation  so  that 


j,(K)  . n(K)^(R) 


(4-19) 


where  j**''  is  also  a row  matrix.  Multiplying  Eg.  (4-18) 
from  the  left  by  integrating  over  electronic  variables 

and  using  Eg.  (4-19)  gives 


(4-20)  becomes 


h<"’  . 


,(K) 


so  that  on  comparing  Egs.  (4-17)  and  (4-22)  one 


(4-23) 


In  an  analogous  fashion  one  arrives  at 


jOCL) 


Using  Egs.  (4-13)  s 


(4-7)  anii  (4-4)  results 


4r 


re 


(4-26) 


The  important  results  of  this  formalism  are  given  by  Egs. 
(4-261  and  14-27).  The  adiabatic  energies  are  found  in 
part  by  solving  the  eigenvalue  problem  for  the  atoms  and 
diatofuicB  that  malce  up  the  molecule.  If  the  eigenvalues  and 
eigenvectors  for  the  fragments  are  )cnown,the  only  electronic 
calculations  that  may  have  to  be  done  are  the  overlap 
integrals. 
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4*5  The  Eigenvalues 

This  section  shall  develop  the  DIM  procedure  to  the 
molecular  ion.  Since  the  results  are  already  available 
in  the  literature  (Ellison,  Huff  and  Patel,  1953)  many  of  the 
details  will  be  omitted.  The  approach  to  be  used  for  this 
system  will  be  itore  physical  in  nature  that  the  general 
outline  presented  in  the  previous  section.  This  will  hope- 
fully clarify  soma  of  the  correlations  between  the  DIM  and 
valence  bond  procedures. 

The  Hj  system  consists  of  three  protons  which  will 
be  labeled  A,  B and  C and  two  electrons.  An  electronic 
function  centered  on  nucleus  A with  an  alpha  spin  will  be 
designated  as  "a"  and  for  beta  spin  "a".  Electronic 
functions  centered  on  nuclei  B and  C are  designated  as  b(E) 
and  c(c)  respectively.  Instead  of  using  a simple  product 
of  atomic  orbitals  for  the  electronic  basis  it  Is  more  con- 
venient to  use  canonical  valence  bond  structures.  There 
are  three  such  structures  for  H,"''  and  they  are  given  by 

*)  = |ah|-|ab|  , (4-28) 

- |be|-|bc|  (4-29) 

and 


(4-30) 
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where  4:  end  63  represent  bonds  between  AE,  BC,  and  AC 
respectively  and  [ — 1 represents  a determinant.  Since  the 
electrons  are  assigned  differently  for  each  structure,  the 
Hamiltonian  will  have  a different  form  for  each  structure. 
The  structures  4i,  4i  and  43  correspond  to  the  Hamiltonians 

H - Hj^g  t ■ “a  ■ "b  ' (4-31) 

H - iigc  + + Hgj.  - ilg  - H|-  (4-32) 


where  AB,  BC  and  AC  refer  to  hydrogen  snolecules,  A,  B and  c 
to  hydrogen  atoms  and  ab,  be  and  ac  to  hydrogen  molecular 
ions.  One  notices  from  Egs.  14-31)  through  (4-33)  and  from 
the  discussion  in  ths  previous  section  that  one  needs  valence 
bond  functions  that  are  eigenfunctions  of  the  fragment 
Hamiltonians.  For  the  isolated  diatomic  AB  one  has 

4“  = |aS  |-  |ao  I , (4-34) 


(a  + b) 


(4-36) 


(4-37) 


where  the  notation  for  the  Hamiltonian  alao  applies  here. 
The  first  function,  Eq.  (4-34),  corresponds  to  the  H,  bond 


where  is  the  ground  state  energy,  ‘E*,  for  H-  which 
depends  on  the  distance  between  nuclei  A and  B.  The  func- 
tions defined  by  Egs.  (4-35)  and  (4-36)  correspond  to  a 


where  i is  either  1 or  2 and  Ei  is  the  ground  state 
energy,  of  which  is  of  course  dependent  on  the 

distance  between  nuclei  A and  B.  The  final  two  functions 
given  by  Egs.  14-37)  and  (4-36)  correspond  to  the  first 
excited  state  of  so  that 


where  i is  either  3 


corresponds 
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excited  state,  of  as  a function  of  the  distance 

between  A and  B.  One  notices  that  Eqs.  (4-35)  through 
(4-30)  furnish  a ineans  of  expressing  the  atomic  orbitals  in 
terms  of  the  valence  )>ond  functions.  For  example,  from 
Eqs.  (4-35)  and  (4-37)  one  has  that 

a . {*?*’  * »f‘>)/2  . (,-42) 

Similar  results  are  obtained  for  the  isolated  diatomlcs 
BC  and  AC.  As  an  example  of  how  to  construct  matrix  ele- 
ments consider 


AjHag^i  * AjH^^ (a5-ab)  , (4-43) 

From  an  analogous  procedure  that  led  to  Eq.  (4-42)  one 
obtains 

AzH^gJi  • AzH^^I(«i°  )B  - (4j®  4 4?°)b]/2  = 

* ae  ac 

Ajt(Ei  »i  4 - (Ei'*f®  4 Ei°»?°)b]/2  . 

Re-expressing  the  valence  bond  states  in  terms  of  atomic 
orbitals,  rearranging  the  terms  and  carrying  through  the 
antisymmetrisation  operation  results  in 

AjHgc*‘  * + «il  * C?®(4,  - tj)J/2  . (4-45, 


Carrying  through 


algebra  for  the  other  cases  results 


AsH^*i  = * «i)  + Ej°(*i  - «s)3/2  , 

* bj)  + e“{«j  - tiU/2  , 

+ «s)  + - »tl]/2 

AsHbc**  “ IS?®(«i  t *j)  + E?°(*s  - *1)1/2  , 

♦ 4,)  + Bf**!*,  - *8)1/2  . 

* E^*1  » 

A2Hjp»j  - 


(4-46) 

(4-47) 

(4-49) 

(4-49) 

(4-50) 

(4-Sl) 

(4-52) 


AjEjU;*!  = E?*^»J  . (4-53) 

Combining  Eqs.  (4-45)  through  (4-53)  with  Eqs.  (4-31) 
through  (4-33)  results  in  the  matrix  elements  of  the 
Haiultonlan  being  defined  by 


<4.|Hi*i>  = [(2E?®  ♦ E?'  f Ef= 


4Eg)S,. 


(Ef°  - Ef®)S|; 


(E,  - El  )Sul/2 


(4-54) 
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= [(Et°  - e?“)S,i  -I-  (2Ef‘=  -f  E?®  * Ef®  + sf’ 

•*  - <SglS,2  t (E?^  - S?^)Si3l/2  , (4-55) 

= HE?®  - E?®)Si,  t (E?*’  - eP>S,2 

♦ (2Et^^  ♦ E?®  + E^®  + E?^  + E?^  - 4Eg)S,jl/2  , (4-56) 

<*2|H|4,>  • K2B?®  + E?°  -►  B?®  + E?®  ♦ B?®  - 4Ej,)Sj, 

♦ (E?®  - Ef®)Sj2  + (B?®  - Et®)SijJ/2  , 14-57) 

<»2|H|*8>  - !(E?®  - E?®)S,i  -f  I2E?'^  + E?®  + E?®  + E?” 

7 £?“  - 4Eg)S2j  4-  (E?*’  - B?*^)S,3J/2  , (4-58) 

<*j|H|»3>  - KB?®  - E?®)Sji  4-  |E?^  - E?^)Sjj 

♦ llE?'^  + E?®  + E?®  ♦ E?*’  •*  - 4Eg)Sj3]/2  , (4-59) 

. (I2E?®  * B?®  4-  E?®  + E?®  t E?®  - 4B^)S,, 

♦ (E?®  - E?®)Ss2  + (E?®  - E?®)Sj,l/2  . (4-60) 

<*.|H|»,>  = KB?®  - E?®)S,3  + I2E?‘^  -I-  B?®  ♦ Ef® 

♦ E?*’  ♦ fif*  - 4Ejj)Sj2  -f  (E?*’  - E?*’)S,il/2  (4-61) 


<*i|H|*3>  = (IS?®  - E?°)Sj,  + (£?*’  - 

•<•  t •*•  E?®  ♦ S?*  + - 4E(j)Sj,l/2  (4-621 


(4-63) 


One  notices  at  tliis  point  that  the  main  objective  of  DIM 
has  been  obtained,  i.e.  the  Hamiltonian  is  expressed  in  ten 
of  isolated  fragment  energies  and  overlaps.  The  overlaps 
axe  easily  found  to  be  given  by 

Sti  - 2(1  + <a|b>M  , (4-641 

Sj!  - 2(1  t <b[c>‘)  , (4-651 

Sjj  - 2(1  t <a|c>’)  , (4-661 


Sis  - Ssi  - 2(<a|c>  + <a|bxb|c>)  , (4-67) 

Si!  - S,i  - 2(<b|o>  + <a|exa|b>)  , (4-68) 


Si!  = 2<<a|b> 


i|  cxbj  c>) 


(4-69) 
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Host  applications  of  DIM  have  asBuned  that  atomic  functions 
centered  on  different  nuclei  are  orthogonal.  This  will 
also  be  assusved  here  so  that  the  eleisents  for  the  overlap 
matrix  are  given  by 


’ij 


With  the  overlap  defined  by  Eg.  C4-70),  the  matrix  elements 
of  the  Kainiltonian  oonslderably  simplify  and  one  obtains 


Hii 


E?°  + Es®  - , (4-71) 


(4-73) 


(4-74) 


Hi  I = H j 1 


E?®  - Ei“ 


Kii 


(4-76) 


J-77) 


One  should  mention  thet  the  KamiltoniaQ  with  ntatrix  elements 
defined  by  Bgs.  (4-54)  through  (4-63)  is  not  Hemitian. 

Howevet  with  the  oveclap  defined  by  Eg.  (4-70)  the  resulting 
Hamiltonian  becomes  Kermitian.  In  other  words,  the  assumption 
for  the  overlap  eliminates  the  non-Hermitian  terms  in  Che 
Hamiltonian. 

Expressing  the  eigenfunctions  of  Che  total  Hamiltonian 
as  linear  combinations  of  the  canonical  valence  bond  function 


(4-70) 


one  obtains  in  accordance  with  Eg.  14-4)  the  following 
matrix  equation 


Hr  - S[E  (4-79) 

where  the  elements  of  Che  diagonal  matrix  E correspond  to 
the  adiabatic  eigenvalues  and  Che  matrix  elements  of  K were 
defined  in  Eg.  (4-77).  As  la  well  known,  the  eigenvalues 
are  found  by  solving 

[ (H  - ES) I a 0 (4-00) 

where  again  | [ indicates  a determinant.  From  Eg.  (4-80) 

it  is  straightforward  to  show  that  the  eigenvalues  ace 


given  by 


- HazHai  - BiiHaa  - HiiHsa)^  + H i iHaaRa a+2Ri aHa 3H1 a 


- HiiHsj 


HlaHas 


vhich  is  a third  order  polynomial  whoae  roote  furnish  the 
three  adiabatic  energies.  It  is  a straightforward  exercise 
to  determine  these  roots  analytically  (Selby,  1964) . The 
lowest  two  roots  correspond  to  the  two  lowest  surfaces  for 
Hj  while  the  other  root  corresponds  to  an  excited  surface. 
This  excited  surface  is  energetically  inaccessible  and  does 
not  couple  to  the  lower  surfaces.  Because  of  this,  the 
excited  surface  will  not  be  considered  so  that  the  problem 
involves  only  two  surfaces. 

4-4  The  Eioenfunctions 

£ven  though  the  adiabatic  energies  were  found  without 
determining  the  eigenvectors,  these  eigenvectors  are 
necessary  to  form  other  quantities  such  as  non-adiabatlo 
couplings.  As  was  mentioned  in  the  previous  section  the 
upper  excited  surface  does  not  couple  to  the  two  lower 
surfaces  so  that  one  has  in  essence  a two-surface  system. 

In  the  following,  the  eigenvalue  and  eigenvector  of  the 
lowest  surface  will  be  designated  as  Ei  and  Ti  respectively. 
This  is  the  surface  that  asymptotically  corresponds  to  H.,. 
The  other  surface  which  asymptotically  eorredponds  to  B* 
will  have  its  eigenvalue  and  eigenvector  designated  as 


El  and  respectively.  Since  fi  and  I's  are  linear  combl- 
nations  of  the  three  canonical  valence  bond  functions, 
there  are  a total  of  six  coefficients  to  be  deCermined. 

From  Eq.  (4-79)  one  has  that 

(Hii  - 2E]]rii  + HisTji  + HiiTj!  * 0 , (4-82) 


Hiir,,  + (Hjj  - 2E))r2i  + HsaTj,  = 0 , (4-83) 


Hsifii  + Hi2rji  + (Hj]  - 2Ei)ri,  - 0 . (4-84) 


Defininq 


Bi  = h5.  - (Hii  - 2Ei)(Hi,  - 2Ei)  , (4-05) 

Si  “ Hij(H22  - 2Ei)  - Hi2Hjj  , (4-86) 


(4-87) 


it  is  not  difficult  to  show  that 


= Si 


(4-89) 


found  by  requiring 


normalized. 


The  coefficients  of  »,  are  determined  by  solving 

[Hii  - 2Ej)r,j  -<•  H,,ru  + = O , (4-91) 

HuTij  + (Hjs  - 2E,)rzj  + = 0 , (4-921 


ri,r,j  + rj,r,s  + r,,r„  = o (4-93) 

where  Eq.  (4-93)  insures  that  the  eigenfunctions  are 
orthogonal . Defining 


(4-94) 


t-95) 


TsiHia  - TiiHt 3 


63  = 


- 2Ej)  - 


r,3  - B,r,2 


l4-9«) 


r,i  • BjCaTii 


(4-97) 


The  value  o£  Tib  la  obtained  from  the  normaliaation  of 


One  notices  that  the  dependence  of  these  coefficients  on 
nuclear  distance  is  contained  solely  in  the  adiabatic 
eigenvalues.  Since  these  are  given  analytically  by 
Eq.  (4-81)  these  eigenvalues  are  also  analytical  functions 
of  the  internuclear  variables.  As  will  be  seen  in  the  next 
section,  this  will  considerably  simplify  the  problein  of 
determining  the  non-adiabatic  couplings. 

4-5  Non-adiabatic  Couplings 

As  was  mentioned  previously,  one  of  the  advantages  of 
DIM  is  that  it  is  relatively  easy  to  obtain  the 


fij  • 


non-adlabatic  couplings.  Since  eigenfunctions  obtained  from 
an  ab  initio  or  other  semi-empirical  approach  are  not 
analytic  functions  of  the  nuclear  variablee.  at  least  two 
complete  calculations  would  have  to  be  done  to  obtain  the 
coupling  of  one  configuration  of  the  nuclei,  in  DIM. 
however,  the  eigenfunctions  depend  analytically  on  the 
nuclear  variables  so  that  the  couplings  are  readily  ob- 
tained functions  of  the  internuclear  variables. 

to  simplify  matters,  the  momentum  operator  will  be 
assumed  to  have  the  form  of  Eg.  (2-4) . Purthermore . only 
one  component  of  will  be  considered  which  will  be 
designated  as  . since  from  Bg.  (2-19)  the  couplings  from 
can  be  expressed  in  terns  of  the  couplings  from  3/3X, 
this  section  shall  only  treat  the  couplings  that  cone  from 
^ . The  matrix  elements  of  interest  are  given  by 

(4-99) 

where  9^  was  defined  in  Eg.  (4-70).  One  can  easily  show  by 
tajting  the  partial  derivative  of  with  respect  to  X 

that  the  matrix  of  these  couplings  is  antisymmetric  with 
diagonal  elements  equal  to  zero. 

Two  methods  for  obtaining  the  couplings  given  in 
Eg.  (4-99)  will  be  presented  here.  The  first  method  will 
obtain  the  couplings  directly  from  the  expansion  coeffi- 
oients  which  depend  analytically  on  X.  The  second  method 
will  employ  partial  derivatives  of  the  Hamiltonian  matrix 


with  respect  to  X and  It  will  be  seen  that  this  may  have 
some  advantages . 

In  general,  the  coupling  matrix  whose  elements  are 


where  its  antisymmetric  property  has  been  used.  From 
Eg.  (4-78)  and  with  the  overlap  given  by  Eg.  (4-70)  one 


+ rji|jjr!j)  . (4-ioi) 


This  may  at  first  glance  seen  to  be  a rather  simple  ex- 
pression. However  one  notices  from  the  equations  for  the 
expansion  coefficients  given  in  the  previous  section  that 
they  are  fairly  complicated  expressions  in  terms  of  the 
matrix  elements  of  the  Hamiltonian  and  eigenvalues.  Analy- 
tical expressions  for  the  partial  derivatives  of  these 
coefficients  may  be  very  complicated  so  that  numerical 
procedures  may  be  more  appropriate  for  calculating  Eg. 


defined  by  Eg.  (4-99)  is  given  by 


(4-100) 
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The  other  method  starts  from  the  time  independent 
Schrddinger  equation.  Tahing  the  partial  derivative  with 
respect  to  X of  this  equation  results  in 

+ « lx  ' 'lx  lx  • (<-102) 


Multiplying  from  the  left  by  and  integrating  over 
electronic  coordinates  results  in 


<»i llxifj’ 


^f.Klx  h1|y»> 


(4-103) 


where  the  Hernitian  property  of  H has  been  used.  Using 
Eq.  (4-78)  in  Bq.  (4-103)  results  in 

<til|xl»s>  = (r,,(r,j|j(  H,,  + r,,|3f  h,j  + a,,) 

+ Til  !ru|j  Hj,  ♦ r,,|3j  B,,  + Tssf^j  H,,) 

* r„(r,j|j  H,i  + rjs|jf  a„  + h„)}  (e,  - s,)"^. 

(4-104) 

Equation  (4-104)  may  seem  to  be  more  complicated  than  Eq. 
(4-101)  but  one  notices  that  the  matrix  elements  of  the  Hamil- 
tonian are  linear  functions  of  the  fragment  energies.  If  these 
energies  are  given  by  rather  simple  expressions,  Eq.  (4-104) 
may  be  determined  analytically.  This  could  bypass  the  use 


of  numecical  methods  to  determine  the  couplings  or  could 
provide  a means  to  check  their  accuracy. 

All  the  essential  ingredients  from  the  electronic 
problem  to  do  a scattering  calculation  have  now  been  ob- 
tained. The  adiabatic  surfaces  are  obtained  by  finding  the 
lowest  two  roots  of  Eg.  (4-81).  The  couplings  are  given  by 
Eqs.  (4-1011  or  (4-104)  with  the  expansion  coefficients 
obtained  from  Egs.  (4-S8)  through  (4-90)  and  Egs.  (4-96) 
through  (4-98).  The  numerical  results  for  the  surfaces  and 
couplings  will  be  presented  in  chapter  six. 


CHAPTER  5 

THE  NOCI^AR  PROBLEM  FOR  COLtINEAR  H,* 

5-1  Hvperapherical  Coordinates 

Hypsrspherical  coordinates  were  originally  introduced 
by  Delves  (Delves,  1950).  Recently  they  have  found  in- 
creasing use  in  guantuiD  mechanical  calculations  that 
Involve  rearrangement  processes  (Smith,  1962;  Kuppentiann, 
1975i  Hauke,  Hanz  and  Rdmelt,  19801.  The  reason  for  this 
is  that  they  describe  both  reactants  and  products  equally 
well.  Using  these  coordinates  avoids  the  necessity  of 
having  to  make  changes  from,  for  example,  one  set  of  Jacobi 
coordinates  that  describe  the  Initial  state  to  another  set 
that  describe  the  final  state  during  the  collision.  In- 
stead, the  problem  is  solved  by  matching,  in  the  asymptotic 
regions,  the  solutions  in  terms  of  hyperspherical  coordi- 
nates to  the  solutions  of  the  asymptotic  channel  Hamilto- 
nians. Thus  the  collision  event  is  described  by  one  set 
of  coordinates. 

This  section  is  concerned  with  developing  the  quantum 
mechanical  Hamiltonian  operator  in  hyperspherical  coordi- 
nates for  a collinear  system  of  three  nuclei.  Although 
the  equations  of  motion  developed  in  chapter  three  do  not 
require  the  formation  of  this 


operator,  constructing 
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helpful  for  obtaining  a more  ccnplete  picture  of  the 
scatterin9  problem. 

To  begin  with,  the  positions  of  the  nuclei  are 
designated  by  z,,  and  z„  defined  in  a center  of  mass 
coordinate  system  such  that  z , -«  < z < is  the  position 
of  nucleus  o relative  to  the  center  of  mass.  In  the 
collinear  case  there  are  two  arrangement  channels  given 
by  <A,BC)  and  (AB,C)  which  will  be  designated  as  channels 
1 and  3 respectively.  The  Jacobi  coordinates  are  defined 
for  channel  1 as 


and  for  channel  3 as 


where  m^^,  and  m^.  are  the  masses  of  nuclei  A,  B and  C 
respectively.  The  Hamiltonian  in  channel  1 is  given  as 
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■°a‘°b  * "‘c’ 

A,BC  i»^  -mig  4- 


[5-7) 


A more  symmetric  form  for  the  kinetic  energy  operator  can 
be  obtained  by  introducing  mass  weighted  Jacobi  coordinates 
defined  in  channel  1 as 


"a,bc, 

"^c 


where  the  channel  subscript  has  been  suppressed.  In  these 
variables  the  Hamiltonian  becomes 

“ * 5^  ‘dH-  * 5yT>  + 

Motivation  for  introducing  hyperapherical  coordinates  is 
provided  by  considering  the  eigenvalue  problem  of  the 


IIS 


kinetic  energy  operation  in  Eg.  (5-10).  The  solutions  of 
this  problem  can  be  written  as 


and  E is  the  continuous  eigenvalue.  Eguation  (5-11)  points 
out  one  of  the  fundamental  differences  between  the  two  and 
three  body  problems.  Since  ki  and  kj  are  continuous  vari- 
ables. the  three  body  problem  in  this  set  of  coordinates 
requires  two  sets  of  functions  that  depend  on  continuous 
variables.  The  two  body  problem,  however,  only  requires 
one  such  set  of  functions.  The  following  change  of  vari- 
ables will  reduce  this  to  one  set  for  the  three  body  problem 
thus  making  it  formally  equivalent  to  a two  body  problem. 

Hyperspherical  coordinates  for  the  case  of  three 
collinear  nuclei  are  defined  as 


ikiX  ikjY 


(5-11) 


(5-12) 


(5-13) 


(5-14) 
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The  Hamiltonien  in  this  set  of  coordinates  becomes 

“ “^4  Ir''  lr>  ♦ ^ If’ 

Considering  the  eigenvalue  problem  for  the  kinetic  energy 
operator  one  has 


(5-17) 


Expressing  the  solution  as  a product  x (r)  S(*)  of  a function 
of  r and  a function  of  e results  in 

!?<'■  Ir’ ^ If  e(»i  • <5-19] 


sides  of  Eq. (5-191  equal 


(5-20) 


r If  ' It 


1^  S(»)  = - Ae(t)  . (5-21) 

The  solutions  of  Eq.  15-21)  can  be  written  as  a linear 
combination  of  sine  and  cosine  functions.  Requiting  that 
0(0)  ■ ® gives 

*n  * (5-22) 


0„(t)  - sin  (5-23) 

where  is  a constant  that  could  be  determined  through 
normalization.  The  full  solution  is  expressed  as 

♦ 1 *„(r)  0jj(O)  (5-24) 

where  satisfies  Eq.  (5-20)  with  A equal  to  A . Equation 
(5-24)  illustrates  one  of  the  advantages  of  using  these 
coordinates  in  that  the  ^ dependence  is  given  by  a denumerable 
(i.e.  countable)  set  of  functions.  This  formally  reduces  the 
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problem  to  a two  body  problem  as  is  apparent  by  noticing 
the  similarity  of  Eq.  (5-201  to  the  radial  part  of  the  two 
body  problem  in  the  center  of  mass  ccordinate  system. 

It  is  convenient  to  further  define  a variable  0 such 


Its  relation  to  f is  given  by 


so  that  the  Hamiltonian  becomes 

u . fi*  ,1  a ,3  ^ ¥*  a*  ^ ^ S , a,  ,e  a,, 

" <r  ?r  5?  * * ‘’al'’"'®’  • ‘^-2” 

One  of  the  reasons  for  introducing  this  angle  is  that  the 
regions  of  0 that  correspond  to  breakup  and  rearrangement 
do  not  depend  on  the  masses  of  the  nuclei. 

Another  advantage  of  using  these  variables  can  be  seen 
by  considering  Bqs.  (5-5)  and  (5-27).  The  Hamiltonian 
given  by  Eg.  (5-5)  is  channel  dependent.  Its  coordinates 
are  particularly  useful  for  describing  the  asymptotic  states 
in  channel  1.  If  a rearrangement  occurs,  it  becomes 
necessary  to  transform  to  another  set  of  coordinates  that 
are  useful  for  describing  the  final  states.  This  introduces 
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practical  problems  such  as  when  to  make  the  transformation. 
Also  since  eigenstates  corresponding  to  different  channels 
are  not  necessarily  orthogonal,  overcompleteness  can  occur. 

These  difficulties  are  avoided  by  using  hyperspherical 
coordinates.  The  collision  process  is  handled  by  matching 
the  eigenstates  of  the  Hamiltonian  in  terms  of  hyperspherical 
coordinates  to  those  of  the  initial  channel  Hamiltonian  in 
the  Jacobi  coordinates  of  the  asymptotic  region.  Hyper- 
spherical coordinates  are  used  throughout  the  oollision. 

When  rearrangements  occur,  0 simply  goes  from  small  values 
for  A + BC  to  values  that  ate  close  to  i for  AB  + C in  the 
asymptotic  region.  Breakups  would  correspond  to  large  values 
of  r,  and  9 In  the  region  of  s/2.  Thus  all  possible  products 
are  described  by  one  set  of  coordinates  so  that  it  is  no 
longer  necessary  to  change  coordinate  systems  during  the 
collision  process.  Also,  since  the  eigenstates  that  are 
functions  of  9 are  orthogonal  and  complete,  overcompleteness 

5-2  The  Hamiltonian  in  the  Almost  Adiabatic  Hepresentation 
This  section  shall  be  concerned  with  obtaining  an 
expression  for  the  Hamiltonian  defined  by  Eg.  (3-59)  in 
hyperspherical  coordinates  with  the  electronic  basis  satis- 
fying Eq.  (2-221.  The  momentum  in  terras  of  hyperspherical 
coordinates  Is  given  by 


(5-291 
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where  is  a unit  vector  in  the  direction  of  r,  Ug  is  a 
unit  vector  perpendicular  to  u and 


(5-29) 


e coordinates, 


e Hasiiltonian  becoraes 


” ' Tn  * amr'ici  * ” (5-30) 

where  for  brevity  the  subsoript  BC  for  the  masses  has  been 
suppressed. 

Expressions  for  the  LapZacian,  divergence  and  gradient 
in  this  coordinate  system  will  be  needed  later  on.  in 


_ 3 


r f « (5-31) 


where  f is  an  arbitrary  function  of  r and  6 and  P is  an 
arbitrary  vector  with 


components 


E*  & i * h’  If  f * iAr  s*  |jy ; * '|j  p 

• F^'li  £’|  'K  =1  * I fii  S • E’lfe  gi; 

• FiW  S* 'll  Si'S  ■•  '3- 


m;Ii,  - <M5;ltj> 


(5-371 
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f COS  »(r,6)  sin  l(r,e)| 

* [-8inX(r,91  cos  »(r,6)J 

where  the  electronic  basis  is  assumed  to  be  real.  By  using 
Egs.  (5-36)  and  (5-38)  in  Eq.  (5-34),  it  is  not  difficult 
to  show  that  the  transformation  angle  satisfies 

3*1  1 31  1 3*1  “^r  *'^r  1 

3?T  ? Jr  * JP’  * T * ~ * 7^  J3“  ° • <5*39) 

Equation  (5-35)  is  a second  order  partial  differential  equa 
tion  which  could  in  principle  be  numerically  solved.  This 
however  far  from  being  a trivial  tas)c.  It  becomes  more 
convenient  to  mahe  a transformation  into  time.  This  is 
accomplished  by  defining 

l(r(t),9(t))  - r(t)  . (5-40) 


By  taking  the  first  and  second  total  derivatives  of  T with 


dr(t) 


(5-42) 


Dsin9  Eqs.  (5-41)  through  (5-44)  in  Bq.  (S-39)  leads,  with 
some  rearrangement,  to 


S{r,8) 


Defining 


dA(t) 


A(t) 


|T{t> 


A(t)  - 0 


with  boundary  conditions 


(5-51) 


where  is  the  initial  time.  Bga.  (5-46)  and  (S-49)  are 
first  order  differential  equations.  They  can  be  solved 
along  with  Hamiltonian's  equation  for  the  nuclear  variables 
to  give  the  transformation  angle  as  a function  of  time. 

5-2b  The  Equations  of  Motion 

In  order  to  determine  the  equations  of  motion  an 
expression  for  the  Hamiltonian  must  be  obtained.  Prom 
Bq.  (5-30)  one  sees  that  the  Hamiltonian  is  expressed  as 
the  sum  of  the  kinetic  energy  and  a potential  energy  (7. 
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The  expression  for  N depends  on  the  electronic  representation 
and  on  the  form  chosen  for  the  nuclear  expansion  coefficients. 
Defining 


there  are  three  terms  that  need  to  be  evaluated.  The 
potential  is  given  as 


where  is  the  adiabatic  potential  matrix.  ?or  the  2x2  case 

Vii  - cos‘r(t)  + V^2  sin*r(t)  (5-S4) 

’'22  * ''11  si>'’''<t)  + V®2  cos*r(t)  (5-55) 


'^12  “ ''21  ■ •''11  * ''22’  cosr(t)sinr{t)  (5-56) 

where  and  are  defined  In  Eq.  (2-11)  and  correspond 
to  the  two  lowest  surfaces  of  H.  . 

Using  Eqs.  (2-21)  and  (5-33)  in  Eq.  (3-23)  leads  to 
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and  was  defined  by  Bq.  (5-36).  with  C defined  by 
(5-30),  one  obtains  after  transforroing  into  time 


3t  J 


(5-59) 


Prom  Egs.  (5-20)  and  (5-59)  the  second  term  in  the  expression 


(5-61) 


where  Eg.  (5-49)  has  been  used.  In  a similar  manner  one  obtains 


,A(t) 


5-62) 
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using  Eqs.  {5-53),  (5-61)  and  (5-62)  in  Eg.  (5-52)  and 
substituting  the  result  into  Eg.  (5-301,  one  obtains  an 
expression  for  the  Hamiltonian  of  the  form 


- ’’r  ^6  ra  '^9* 

“ “ Ira  * 2itir!r>  * r’‘l~  * * p|-  ''’<«<  * EFT?? 

2d*  j s 

* A(t)  ♦ ^ A(t)«)  4 (c;c„)-lc;(V  4 P^v* 

i>SN 

where  the  subscript  B has  been  added  to  the  nuclear  expansion 
coefficients  to  avoid  confusion  with  the  transformation 

The  expression  for  ff  given  in  Eg.  (5-64)  is  quite 
complicated  and  the  final  expression  depends  on  the  form 
chosen  for  the  nuclear  expansion  coefficients.  It  is 
straightforward,  given  the  form  of  the  nuclear  expansion 
coefficients,  to  determine  the  final  font  of  the  Hamiltonian. 
However,  since  the  calculations  were  performed  in  a different 
electronic  basis,  they  will  not  be  presented  here.  If  this 


representation  were  used,  one  would  need  to  simultaneously 
solve  a set  of  ten  first  order  differential  equations. 

Pour  of  them  would  correspond  to  the  nuclear  positions  and 
itomenta  and  four  would  determine  the  nuclear  expansion 
coefficients.  The  other  two  given  by  Eq.  (5-48)  and  (5-49) 
are  needed  to  determine  the  transformation  angle. 

5-3  The  Hamiltonian  in  the  Diabetic  Representation 

The  diabatic  representation  is  defined  by  requiring  that 
the  electronic  basis  satisfy  Eg.  (2-13).  There  ace  two  main 
advantages  in  using  this  representation.  The  first  is  that 
the  transformation  angle  is  much  easier  to  determine.  The 
second  is  that  this  representation  leads  to  a rather  simple 
form  for  the  Hamiltonian.  Another  apparent  advantage  is 
that  it  seems  to  help,  from  a computational  point  of  view, 
if  the  couplings  do  not  appear  explicitly.  The  rest  of  this 
chapter  will  develop  the  electronic  transformation  and 
Hamiltonian  in  this  representation. 

S-3a  The  Electronic  Transformation 

As  was  discussed  in  section  (2-2b),  the  electronic 
transformation  satisfies  Eg.  (2-15).  In  hyperspherical 
coordinates,  one  has  that 


where  the  transformation  was 


chosen  of 


sin(X)) 


and  was  defined  in  Eq.  (5-37).  By  comparing  the  partial 
derivative  of  d®  with  respect  to  9 to  the  partial  derivative 
of  dg  with  respect  to  r,  one  finds  that 


(5-69) 


ii(r,9) 


(5-70) 


5-67)  in  (5-70)  results  ii 


= - a*(r,8) 


a»,r.v  = 


>lr,ei  . MvV  - (‘ 


f Eq.  (5-73)  BO 


,=*=,-‘c*.«= 


(5-7S) 


(S-77) 


where  X is  given  by  Bq.  (5-73) . II 
the  nuclear  expansion  coefficients 


the  synooetric  form  for 


(S-78) 


(5-79) 


This  essentially  completes  the  development  in  that  the 
Bamiltonian  is  found  by  using  Bq.  (5-80)  in  Eq.  (S-30). 

Doing  this  results  in  a set  of  eight  first  order  differential 
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e^ations  in  time  that  need  to  be  solved  simultaneously. 
Explicitly,  these  equations  are 


as 

ar 


a*! 

3t“ 


3H l,„d  „ . „d  „ , 

ax^  ■ K'^ll'^1  ''12  2' 


■“2  _ 3H_ 
3t  3^2 


(5-ai) 


(5-82) 


(5-63) 


(S-84) 


(5-35) 


(5-67) 


(5-88) 


The  analytical  expressions  that  determine  P and  hai. 
been  omitted  because  of  their  complexity.  In  practice 
is  found  that  they  are  easier  to  determine  numerically. 


CHAPTER  6 

CALCOIATIONS  FOB  H3*' 

6-1  Electronic  Reaults 

This  section  will  be  concerned  with  the  numerical 
aspects  of  applying  the  equations  developed  in  chapter  four 
to  the  collinear  system.  It  was  shown  in  chapter  four 
that  the  method  of  Diatomics  in  Molecules  with  zero  dif- 
ferential overlap  led  to  a solution  of  the  eleotronic 
problem  that  depended  solely  on  tJie  diatomic  potentials. 

The  following  subsection  will  develope  the  diatomic 
potentials  that  are  necessary  for  the  Diatomics  in  .'lolecules 
treatment.  These  will  be  used  in  the  following  subsections 
to  construct  the  potentials  and  couplings  for  the  collinear 
system  which  will  he  given  in  terms  of  hyperspherical 
coordinates. 

6-la  The  Diatonic  Potentials 

AS  was  discussed  in  Section  (2-3),  for  each  pair  of 
nuclei,  three  diatomic  potentials  are  needed  to  detemine 
the  electronic  problem.  These  correspond  to  the  ground 
state  potential,  'z*.  for  Hj  and  the  ground  state,  ‘l*,  and 
first  excited  state,  I*,  potentials  for  Hj*. 

A Morse  function  (Morse,  1929)  was  chosen  to  represent 
the  ground  state  for  Hj.  Explicitly  one  has 
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Fig,  (12)  The  )4orse  ratential,  Bq.  (6-1),  for  the  ground 
state,  of  H2.  The  parameters  are  defined 

in  the  text  following  Eq.  (6-1). 


>E*(r)  - >0(e"  - 2e"  '®'*’  ‘ ‘*’o’)  (6-1) 

where  r is  the  Internuclear  distance.  For  the  well 
depth,  ‘d,  is  4.759  e.V.,  the  position  of  the  minimum,  ‘r  , 
is  1.402  a , and  ‘a  is  1.044  a'^.  The  potential  is  shown 
in  Fig.  (12) . 

The  ground  state  for  H-  was  also  given  the  form  of  a 
Morse  potential  so  that 


where  and  *a^  have  the  values  2.795  e.V.,  2.003  a^, 
and  .7194  a respectively.  An  anti-Morse  potential,  defined 
by 


was  used  to  describe  the  first  excited  state,  *I*,  of  Hj. 
The  values  of  the  parameters  were  chosen  to  be  17.1  e.V., 
.344  a^  and  .6708  a^^  for  *0^,  ^r^  and  *a^  respectively. 
The  values  of  the  parameters  used  in  Egs.  (6-2)  and  (6-3) 
are  essentially  the  same  as  those  used  by  Petersen  and 
Porter  (Petersen  and  Porter,  1967).  These  potentials  are 
shown  in  Fig.  (13) . 

The  use  of  Morse  and  anti-Morse  potentials  may  seem 
to  be  a bit  arbitrary  in  that  more  accurate  diatomic 
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potentials  could  be  easily  obtained.  They  are,  however, 
oonvenient  for  maXing  ooinparisons  in  that  other  theoretical 
studies  (Tully  and  Preston,  1971i  Top  and  Baer,  1977)  have 
used  potentials  of  this  form, 

6-lb  Adiabatic  Potential  Surfaces 

The  adiabatic  potential  energy  surfaces  are  found  by 
solving  for  the  roots  of  the  third  degree  polynomial  defined 
in  Eg.  (4-81) . These  can  be  determined  analytically  (Selby, 
19649  and  the  lowest  two  roots  correspond  to  the  surfaces 
that  are  needed  in  the  scattering  calculations. 

The  diatomic  potentials  that  were  presented  in  the 
previous  section  depend  on  intemuclear  distance.  In  terns 
of  hyperspherical  coordinates,  these  internuclear  distances 
become 


where  r^g  is  the  distance  from  nucleus  a to  nucleus  6 and 
the  mass  terns  were  defined  in  Section  5-1. 
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Non-adiabatic  couplings  due  to  6 


Fig.  (18) 


),  Eq.  16-7)  . 


n-adiabatic  couplings  d 


on<*adiabatic  couplings  due 


(6-7) . 


Fig.  (20) 


Fig.  (21J  Non-adiabatic  couplings 
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Non-adiahatic  couplings  due 


Fig.  (22) 
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Fig.  (23) 


n-adlabdtic  couplings  due 


Bq.  (6-8). 
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Pi?.  {24)  Kon-adiabBtic  couplings  due  to  c.  Xg.  (6-3). 


-adiabatic  couplings 


The  couplings  due  to  6 given  in  Eq.  {G-7)  are  shown  in 
Pigs.  (18)  through  (21)  as  a function  of  8 foe  different 
values  of  r.  Since  the  couplings  are  symmetric  about  v/2r 
only  the  results  for  0<8£t/'2  ace  shown.  For  r equal  to  48^, , 
the  surfaces  are  far  apart  and.  as  one  would  expect,  the 
couplings  are  saiall  except  for  small  values  of  8 . For  small 
values  of  8 the  surfaces  are  highly  repulsive  so  that  these 
oouplings  lie  in  a region  that  is  in  this  study  energetically 
inaccessible.  As  r increases,  the  coupling  becomes  stronger 
and  more  localised  in  the  region  where  the  potentials 
approach  each  other.  For  r equal  to  10  a^,.  the  coupling  is 
large  and  localised  about  the  (>otential  pseudocrossing  and 
eventually,  for  r in  the  asymptotic  region,  the  coupling 
becceies  a delta  function  (Top  and  Baer.  1977]  . 

The  couplings  due  to  r given  by  Eq.  (6-9)  are  shown  in 
Figs.  (22)  through  (26)  for  different  values  of  r.  One 
notices  that  they  have  the  same  general  behavior  as  the 
couplings  due  to  8 but  they  are  also  much  weaker.  This  is 
one  of  the  advantages  in  using  hyperspherical  coordinates 
in  that  since  the  couplings  due  to  r are  much  smaller  than 
those  due  to  8 . it  may  be  possible  to  neglect  them. 

6-Id  model  of  6Ton-adiabatic  Couplings 

The  results  to  be  presented  in  the  next  section  were 
done  in  the  diabetic  electronic  representation.  In  chapter 
five  it  was  shown  that  the  traneformation  from  the  adiabatic 
to  the  diabetic  representation  required  only  the  non- 


adiabatic  couplings  due  to  B . It  was  further  suggested  that 
if  these  couplings  exhibit  a rather  simple  behaviort  it  may 
be  possible  to  model  them  with  simple  functions  so  that  the 
integral  in  Eg.  (S-73)  could  be  determined  analytically. 

In  the  previous  section  it  was  seen  in  Pigs.  (18) 
through  (21)  that  the  couplings  due  to  9 were  fairly  sym- 
metric around  the  pseudocrossing.  A function  with  this 
property  is  the  Gaussian,  so  that  the  couplings  were  chosen 

, - b(r)  (9  - e„(r))^ 

d|(r,6)  = B(r)e  * (6-9) 

where  9^^  is  the  value  of  6 that  corresponds  to  the  maximum 
of  the  coupling.  The  range  of  9 in  Eg.  (6-9)  is  restricted 
to  be  between  rero  and  ir/2.  The  coupling  for  values  of  9 
greater  than  s/2  are  simply  found  by  evaluating  Eg.  (6-9) 
with  9 equal  to  * - 0 . B(r)  was  chosen  to  be 


where  A * .045,  a ■ .SSa^^  and  the  second  term  on  the  right 
hand  side  was  added  to  improve  the  agreement  at  r equal  to 


1S4 

where  a = 1.1684,  ~ -.04215  a^^  and  again  the  last  term 

improved  the  agreement  around  r equal  to  Sa^.  Finally,  for 


= 1.7  and  e - l-lSSa^ 


e derivative  of 


b(r)  i 


. b(r)  I 


that  the  integral  is  still  well  defined.  The  results  of 
the  couplings  defined  by  Eg.  (6-9)  are  compared  to  the 
actual  couplings  in  Figs.  (26)  through  (29).  In  general 
the  agreement  is  not  bad.  The  worst  agreement  occurs  for 
r equal  to  43q  but  the  coupling  is  small  in  this  region 

6-2  Trajectory  Calculations 

Since  the  collision  process  occurs  in  a very  short 
amount  of  time  (10“^*  secs  - lO"^®  sees),  it 
to  choose  a set  of  units  that  reflect  this.  There  ar 
units  that  need  to  be  defined  and  they  correspond  to 
for  mass,  energy,  length  and  time.  A convenient  unit 
length  is  Che  atomic  unit,  a^,  which  is  equal  to  .529 
A natural  choice  for  the  unit  of  mass  is  the  proton  r 
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the  iDodeled  aouplings  { ) 


Pig.  (26) 


Fig.  (27)  A coapariaon  of  the  modeled  oouplinga  ( ) 

due  to  0 to  the  actual  ones  (* ) , r « 6 a . 
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Pig.  A canparieon  of  the  modeled  couplings  i ) 

due  to  0 to  the  actual  ones  ( ) . r = 8 a . 


158 


(88)  A coinpArison  of  tho  modeled  couplings  (---J 

due  to  9 to  the  actual  ones  ( ) , r * 10  a , 
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nass,  n (1.6726  x 10  , and  electron  volts,  e.V. 

Cl. 6022  X 10  ^^ergs) , will  be  chosen  for  the  unit  of  energy. 
The  unit  of  time,  t,  corresponding  to  the  above  units  is 
given  by 


which  gives  t equal  to  5.407  x lO'^^  secs.  The  value  of 
Planck's  constant  in  these  units  is  found  to  be  .121733 

As  will  become  apparent  later,  it  is  necessary  to 
define  transforraations  between  Jacobi  and  hyperspherical 
variables.  In  terms  of  the  hyperspherical  coordinates 
(r,0),  the  Jacobi  coordinates  in  channel  1 (A,BC)  ace 
given  by 


*1  ' („  '•  ) r CQs(t9)  (6-16) 

where  and  were  defined  in  Eqs.  (5-6)  and  (5-7)  and 

r was  defined  by  Eq.  (5-29).  In  tenns  of  Che  Jacobi 
coordinates,  r and  6 are  found  to  be 


(6-17) 


^1,  °^C  )■>) 
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by 


(6-; 


*3  = ^1  - 5^^ 

■ =1^  ■*  ■ vs  -X' ■'  ■ "■”’ 

'-3  ■ S **-b'S  • ’>1  ■ '■!  ' n-m 


to  be  1 and  the  separation  distance  was  set 
1.402  a , The  mass  of  the  inooining  atom  was  1 m and  its 
initial  distance,  Z. , was  chosen  to  be  10  a . With  the 
internal  energy  of  the  diatomic  set  to  zero  and  with  a 
total  energy  of  1 e.v.  the  initial  moraants  ?,  can  be 
determined.  The  momenta  is- of  course  equal  to  zero. 
With  these  initial  conditions  in  terms  of  Jacobi  coordi- 
nates. Eqs.  (6-17),  (6-18),  (6-22)  and  (6-23)  are  used  to 
determine  the  initial  conditions  in  terms  of  hyperspherical 
variables.  The  Hamiltonian  is  found  to  be 

H(r,9)  = 5^4-  cos(re))*2  (6-28) 

where  Eq.  (6-16)  was  used.  The  equations  of  motion  are 
given  by 


(6-31) 
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(6-32) 


Eqa.  (6-29)  through  (6-32)  were  solved  simultaneously  using 
a standard  subroutine  based  on  the  Rung!,  Kutta,  Gill 
procedure  (Shanpine  and  Allen,  1973)  to  do  the  integration. 
The  value  of  n in  these  calculations  was  set  equal  to 


Figs.  (30)  and  (31).  In  Fig.  (30)  the  coordinates  r and  e 
are  plotted  as  a function  of  time.  One  sees  that  r is 
initially  around  11  a^  and  it  steadily  decreases  until  it 
reaches  a turning  point  at  a time  of  about  5.75  t.  After 
the  turning  point,  r increases  as  the  particles  recede. 

The  behavior  of  6 is  also  quite  simple.  Initially  it  is 
small  and  it  increases  to  somewhat  larger  than  2.5  at  the 
turning  point.  After  the  turning  point  it  decreases  bac)c 
to  small  values.  This  behavior  is  what  one  would  expect 
from  a repulsive  potential  in  that  since  rearrangements  and 
complexes  do  not  occur  8 must  initially  and  finally  have 
small  values.  Also  shown  in  Fig.  (19)  is  the  potential  as 
a function  of  tine.  In  the  interaction  tone  it  increases 
to  a value  of  1 e.V.  at  the  turning  point.  This,  of  course, 
coincides  with  the  (cinetic  energy  being  tero. 

In  Fig.  (31)  plots  of  and  versus  time  are  shown. 
Py  initially  is  negative  as  the  particles  approach  each 


The  results  of  these  calculations  are  shown  in 
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fig.  (30) 


t for  the  case  of 
fhe  upper  figure. 


and  e ( ) as  a function  of 

a repulsive  potential,  aZ?'. 
The  potential  is  shown  in'' 


other,  goes  to  aero  at  the  turning  point  and  becomes  positive 
as  the  particles  recede.  P.  is  initially  positive  as  e 
increases,  goes  to  zero  at  the  turning  point  and  becomes 
negative  as  d beccoies  small, 

Although  this  was  a very  simple  test  case,  seme  of  the 
general  behavior  for  these  variables  will  still  be  observed 
for  the  more  complicated  cases,  one  of  the  differences 
between  this  case  and  the  more  complicated  cases  is  the  oscil- 
latory behavior  of  0 and  P,  due  to  the  vibrational  motion 
of  the  diatomio. 


7he  second  test  case  to  be  considered  here  is  the 
collinear  collision  of  a proton  with  a hydrogen  molecule 
with  only  the  lowest  electronic  surface  present.  Since  this 
case  includes  vibrational  motion  and  rearrangements  and 
inelastic  collisions  are  possible,  this  study  should  offer 
some  insight  in  what  to  expect  for  the  behavior  of  the 
hyperspherical  variables  when  two  surfaces  are  present. 

This  system  is  also  convenient  for  determining  the 
guickest  numerical  method  for  integrating  the  equations.  A 
comparison  was  made  between  the  Runge,  Sutta,  Gill  procedure 
used  in  the  previous  test  case  and  a predictor  corrector 
method  based  on  Adam's  formula  (Shampine  and  Gordon,  1975J . 
Trajectories  were  run  with  both  procedures  using  the  same 
initial  conditions  and  error  tolerances.  It  was  found 
that  the  trajectories  obtained  from  the  procedure  based  on 
Adam's  formula  used  about  an  order  of  magnitude  less 
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where  is  the  total  energy,  what  remains  is  to  choose 

the  initial  values  of  the  vibrational  coordinate,  sj,  and 
moToentuiD  , The  value  of  the  initial  vibrational 
momentuia  was  chosen  to  be  either  plus  or  minus  sene  fraction 
of  the  maximum  initial  vibrational  momentum,  which 

is  given  by 


/2mgg(*D  ♦ . !6-35) 

The  initial  values  of  the  coordinate  are  found  to  be 

Sj  - « *r^  - ln{l  + I'l  + a)  (€-36) 


t - 


(6-37) 


and  the  + indicates  initial  values  that  are  less  than  (-) 
or  greater  than  (t)  the  equilibrium  distance. 

The  first  set  of  trajectories  to  be  considered  here 
originated  in  the  n^O  vibrational  level.  From  Eqs.  (6-33) 
and  (6-34)  the  relative  momentum,  P«  , was  found  to  be 
-1,85  m a The  maximum  value  of  the  diatomic  raomentxim. 


was  obtained  from  Eq.  (6-35)  which  gave  0,5189 
The  initial  value  of  the  diatomic  momentum,  p , was 


the  vibrational  coordinate,  of  1.402  a , the  equilibrium 
distance.  The  initial  values  of  the  hyperspherical  variables 
were  found  by  using  Eqs.  (6-17),  (6-18),  (6-22)  and  (6-23). 
The  Hamiltonian  was  given  by  Eq.  (5-30)  with  W equal 
to  the  lowest  electronic  surface  which  was  designated  as 
in  Section  (4-4)  . The  equations  of  motion  were  foxxnd 


to  be 


3t*°  ■ yr  <e-4i) 


where  the  superscript  "a'  indicates  the  adiabatic  surface  and 


Fig. 


f321  Plots  of  r and  as  a function  of  t for  the 

adiabatic  collision  of  with  H2  on  the  lowest 
Bj  surface.  n*0,  -1,9175  e.V.  and 
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Fig.  (3€)  Plots  of  E,  and  P for  the  adiabatic  collision 
of  with  H2  on  t^e  lowest  nt  surface.  n*0. 


175 

where  the  terms  involving  partial  derivatives  were  evaluated 
numerically . 

The  results  of  Eqs.  (6-38)  through  (6-41)  ace  shown  in 
Figs.  (32)  and  (33).  One  notices  tliat  they  exhibit  the  same 
general  behavior  as  in  the  previous  test  case.  The  essential 
difference  is  their  oscillatory  nature  due  to  the  vibra- 
tional motion.  This  is  especially  evident  in  the  trajec- 
tories for  8 and  Pg . One  also  notices  that  initially  6 has 
small  oscillations  but  that  after  the  collision  its 
oscillations  are  noticeably  larger.  This  of  course  corre- 
sponds to  a transfer  of  energy  from  the  translational  to 
vibrational  nodes.  Since  e returns  to  small  values  after 
the  collisionr  this  is  an  example  of  a process  that  doesn't 
involve  rearrangements. 

The  results  are  also  given  in  terms  of  Jacobi  variables 
in  Figs.  (34)  through  (36).  One  sees  in  Fig.  (34)  that 
becomes  constant  both  before  and  after  the  collision.  This 
constant  behavior  served  as  a criteria  for  terminating  the 
calculation.  The  gain  in  vibrational  energy  is  also  evident 
from  the  plots  of  and  given  in  Pig.  (35)  . As  will  be 
seen  shortlyf  the  trajectories  of  given  in  Fig.  (36)  and 
of  serve  as  an  aid  in  understanding  the  collision  process. 

A particularly  informative  quantity  is  the  potential 
energy  as  a function  of  time  shown  in  Fig.  (37).  Before 
the  collision  there  are  small  oscillations  in  due  to  the 
vibrational  motion  of  BC  (for  brevity  the  nuclei  will  be 
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referred  to  as  A,  B and  C and  the  dlatomica  as  AS  and  BC] . 
As  A approaches  there  is  a sharp  decrease  in  followed 
by  a strong  repulsion  when  the  distance  between  A and  B is 


at  a miniouiQ.  As  A and  B separate  there  is  another  large 
decrease  in  and  eventually  A escapes  leaving  BC  with 
greater  vibrational  energy. 

Further  insight  is  obtained  by  considering  Figs.  (35) 
and  (36).  One  sees  in  Pig.  (35)  that  in  the  initial  part 
of  the  collision  around  t - 3.5  i,  BC  is  contracting.  Due 
to  the  large  attraction  between  A and  B,  the  contraction  is 
slower  for  a while.  As  A and  B separate.  BC  further 
contracts  and  then  undergoes  a large  expansion.  One  might 
expect  that  if  BC  were  initially  expanding  in  the  collision, 
the  repulsion  between  B and  C would  be  strong  enough  to 
cause  a rearrangement. 

The  second  set  of  trajectories  shown  in  Figs.  (38) 
through  (44)  were  done  with  the  same  set  of  initial  con- 
ditions that  were  used  tor  the  previous  ones  with  the 
exception  that  the  sign  of  the  diatomic  momentum  was 
changed,  i.e.  P_  * -.5189  mpa.  . Since  0 in  Pig.  (39) 


goee  from  small  values  to  values  close  to  s.  this  is  an 
example  of  a rearrangement.  One  notices  in  Fig.  (42)  that 
is  constant  after  the  collision  and  this  was  used  to 
terminate  the  calculation.  In  practice  and  were 
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Fi9. 


for  the  adlabetic  colll- 
le  lowest  surface, 
and  Pj  - -.S169  mpagT'^. 
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Fig. (42)  Flota  o 


colliaion 
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rig.  (43)  Plots  ot  *3  and  vs.  t for  the  adiabatic 

collision  of  H with  Hy  on  the  lowest  ht  surface. 

''■Of  ®tot“  e.V.  and  = -.5189  tpa„T*^. 
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calculated  throughout  the  collision  and  the  calculation  was 
tive  and  r was  greater  than  some  value  which  was  chosen  to 

In  Fig.  (41)  one  notices  that  during  the  initial  part 
of  the  collision  around  t ■ 3.5  t,  BC  is  expending.  Since 
B is  strongly  attracted  to  A,  this  expansion  is  somewhat 
greater  which  corresponds  to  an  increase  In  the  vibrational 
energy  of  BC.  When  B and  C contract,  they  get  close  enough 
to  undergo  a strong  repulsion  which  gives  rise  to  the  second 
sharp  increase  in  shown  in  Fig.  (44).  Due  to  Che 
increased  vibrational  energy  of  BC,  this  repulsion  is  strong 
enough  so  that  C escapes  and  AB  is  formed. 

The  previous  two  examples  offer  some  insight  into  why 
and  how  rearrangements  cake  place.  The  rearrangement  process 
seems  to  depend  mainly  on  whether  there  is  a substantial 
increase  in  the  vibrational  energy  of  SC  during  the  initial 
part  of  the  collision.  If  BC  is  contracting  during  the 
initial  part  of  Che  collision,  Che  attraction  between  A and 
B causes  this  contraction  to  be  slower  for  a while.  Con- 
sequently there  is  no  initial  increase  in  the  vibrational 
energy  of  BC  So  that  A escapes  and  BC  remains  bound.  If, 
however,  BC  is  initially  expanding  during  the  collision, 
the  attraction  between  A and  B causes  BC  to  expand  to 
greater  distances.  This  enhanced  expansion  corresponds  to 


166 

an  incceasa  in  the  vibrational  energy  of  BC  which  is  great 
enough  for  C to  escape  and  for  AB  to  form. 

Calculations  were  also  done  at  the  same  total  energy 
for  the  n-4  vibrational  level.  These  results  along  with 
the  n-0  results  seem  to  indicate  that  the  mechanism  dis- 
cussed above  is  almost  always  capable  of  predicting  the 
correct  final  nuclear  configuration.  As  is  expected,  it  can 
give  incorrect  predictions  in  oases  where  long  lived 
complexes  are  formed.  It  will  be  seen  later  that,  since 
this  mechanism  only  depends  on  the  initial  part  of  the 
collision,  it  will  still  be  of  use  when  the  upper  surface 
is  included. 

As  is  clearly  seen  in  Fig.  (371,  this  system  has  the 
ability  to  transfer  energy  between  translational  and  vibra- 
tional modes.  Figs.  (35)  and  (43)  help  in  understanding  why 
this  is  so.  One  sees  in  Fig.  (35)  that  as  A is  leaving  the 
collision  region,  around  t = 4.2  t the  diatomic  is  expanding. 
The  attraction  between  A and  B causes  BC  to  expand  to 
greater  distances  which  increases  its  vibrational  energy. 

The  opposite  case  is  shown  in  Fig.  (43)  where,  at  around 
t ■ 4.5  T,  AB  is  contracting.  The  attraction  between  B and 
C slows  down  this  contraction  for  a while  leaving  AB  with  a 
lower  amount  of  vibrational  energy. 

As  a checlc  in  the  calculations,  the  total  energy,  H, 
was  calculated  throughout  the  collision  and  it  was  found 
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to  remain  constant.  Time  reversal  studies  were  also  done 
and  they  led  to  the  same  results. 

6-2b  H3+  Results 

This  section  shall  report  the  results  of  some  typical 
trajectories  that  were  obtained  by  including  both  electronic 
surfaces  for  the  collinear  Hj"*  system.  The  calculations 
were  done  in  the  diabatic  electronic  representation  so  that 
the  results  were  obtained  by  integrating  Eqs.  (5-811 
through  (5-88) . The  electronic  coupling  terms  were  modeled 
according  to  Eg.  (6-9)  so  that  the  angle,  Eg.  (5-73),  that 
determines  the  transformation  between  the  adiabatic  and 
diabatic  representations  is  given  in  terms  of  error  functions 
as  in  Section  (2-5) . 

Common  to  all  trajectories  are  tlie  initial  value  of 
83  and  the  total  energy.  As  in  the  previous  adiabatic  test 
case,  the  initial  value  of  Zi  was  chosen  to  be  14  a and 
the  total  energy  was  set  equal  to  -1.9175  e.v.  (1  e.V.  above 
the  n=4  vibrational  level  of  Hj) . Also  common  to  all 
trajectories  were  the  initial  values  of  the  electronic 
amplitudes.  was  set  equal  to  /Tfi  and  Y^,  Xj  and  Yj  were 
all  chosen  to  be  sero.  From  Eqs.  (3-103),  (3-106)  and 
(3-107)  this  corresponds  to  an  initial  value  of  Pj  « 1 and 
to  initial  values  of  Yj  and  Pj  of  zero.  >3  was  also  set 

The  remaining  initial  conditions  were  chosen  in  a 
fashion  similar  to  the  previous  adiabatic  test  case. 
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Selecting  the  initial  vibrational  energy  level,  the  initial 
internal  energy  ia  found  from  Eq.  (6-331 . This  is  used  in 
Eq.  (6-34)  to  datemine  the  initial  value  of  Choosing 

an  initial  value  of  between  the  turning  points,  the 
initial  value  of  is  given  by 


The  initial  values  of  the  hyperspherical  variables  were 
determined  by  using  Eqs.  (6-17),  (6-18),  (6-22)  and  (6-23). 

The  first  set  of  trajectories  consists  of  an  example 
of  a non-rearrangement  process  that  initiates  in  the  n equal 
to  tero  vibrational  level.  From  Eqs.  (6-33)  and  (6-34) 
the  initial  value  of  P^^  ia  -1.65  mpa^  . Mith  an  initial 

value  of  of  1.45  a^,  Eq.  (6-42)  gives  the  value  of  P^ 
of  .508  nipaji  where  the  plus  sign  was  chosen. 

The  results  of  r,  Pj.  and  8,  Pj  are  shown  in  Figs. 

(45)  and  (46)  respectively.  Their  behavior  is  similar  to 
the  adiabatic  test  case  involving  non-rearrangement. 
Initially  r decreases  monotonically  and  P is  negative  and 
undergoes  small  oscillations.  As  the  particles  collide, 
r reaches  a minimum  and  begins  to  increase  and  P^  goes  from 
negative  to  positive.  For  t greater  than  about  6 t,  r 
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iiiOTiotonically  increases  and  P_  oscillates  around  a positive 
value.  The  final  rate  of  increase  of  r is  less  than  its 
initial  rate  of  decrease  which  signifies  an  increase  in 
vibrational  energy.  This  is  also  seen  by  considering  the 
behavior  of  6 and  P. . Initially  0 is  small  and  undergoes 
small  oscillations.  It  increases  to  about  2 when  the  par- 
ticles are  all  close  together  and  then  returns  to  small 
values  but  with  significantly  larger  oscillations. 

The  results  of  Zj,  and  z^,  are  shown  in  Pigs. 
(47)  and  (48)  respectively.  One  notices  that  for  t greater 
than  6 r>  Pj^  becomes  constant  and  this  was  used  to  terminate 
the  calculation.  One  also  sees  in  Pig.  (48)  that  during  the 
initial  part  of  the  collision,  t oBt,  BC  is  undergoing  a 
contraction.  Since  this  is  a non-rearrangement  process,  the 
meohanism  described  for  the  adiabatic  test  cases  might  apply 
Co  this  system  also.  This  subject  will  be  discussed  again 
a little  later. 

Pig.  (49)  shows  the  time  dependence  of  W.  The  diabatic 
potentials  E^,  E2  and  used  to  calculate  W are  shown  in 
Figs.  (50),  (51)  and  (52)  respectively.  The  adiabatic 
potentials  and  shown  in  Pigs.  (53)  and  (54) 

respectively  and  the  transformation  angle,  1,  is  depicted 
in  Fig.  (55) . One  notices  the  similarity  in  the  figures 
for  W,  and  E^.  Initially  they  undergo  small  oscillations 
characteristic  of  the  lowest  vibrational  level.  As  A 
approaches  B,  there  is  a large  decrease  in  energy  and  the 
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Fig.  (47)  Plots  of  and  vs.  t for  the  non-adiabatlc 
collision  of  H with  on  the  lowest  two 
surfaces,  n-0,  8^^^“ 
and  P^  ”.S08  ™p*(jT“l. 
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Plot  of  H VS.  t for  the  non-sdiabatic 
collision  of  H with  H2  lowest 

surfaces.  n«0,  £^0^—1. 9175  e.V.,  zj«l 
and  P“  ■,506  ra^a^T*^, 


196 


Plot  of  the  diabatic  potential  Bj  va.  t for  the 
collision  of  H*  with  Kj  on  the  lowest  two  Hj 
surfaces.  n^O,  ^rot*  “1.9175  e.V.,  2^=1. 45  a^ 


Pig.  (51) 
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Fig.  (52)  Plot  Of  the  diabatic  potential  vs.  t for  the 
oollialon  of  H*  with  Hj  on  the  lowest  two  H* 
surfaces.  n=0,  6^^^=  -1,9175  e.V.,  tj=1.45a^ 
and  P^  =.509  m a t"1. 
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Fig.  (53)  Plot  of  the  adiabatic  potential  E®  vs.  t for 
the  collision  of  H*  with  H,  on  the  lowest  two 
Hj  surfaces.  n=0,  Etot“  “I- 
and  P^^=.508  »pa^T-i. 
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Pi9.  (55) 


Plot  of  the  txanefonoatlon  angle,  X,  vs 
the  non-ediabatic  oollialon  of  H vitb 


t for 
•1.3175 
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system  is  in  the  well  legion  of  the  lowest  surface.  As  A 
gets  closer  to  B a strong  repulsion  occurs  which  corresponds 
to  Che  sharp  increase  in  the  potentials.  As  A and  B 
separate  there  is  another  large  decrease  in  the  potentials 
as  the  system  goes  through  the  well  region.  The  large  final 
oscillations  Indicate  a significant  increase  in  vibrational 
energy.  On  close  inspection  one  finds  that  initially  and 
finally  Ej  and  E*  coincide  and  that  in  the  collision  region 
is  greater  than  E*  as  it  should  be.  Due  to  the  similarity 
between  E^  and  H after  the  collieion  one  would  expect  a low 
probability  for  electronic  transition. 

The  results  of  and  Pj,  f2  are  shown  in  Figs. 

(56)  and  (57)  respectively  and  fj  and  yj  are  seen  to  oscil- 
late very  rapidly.  The  behavior  of  the  probabilities  can  be 
better  understood  by  considering  the  graphs  of  and 
It  is  well  known  from  the  semi-clasaioal  treatment  of  this 
system  (Tully  and  Preston,  1971)  that  the  derivative  of  the 
amplitude  for  the  upper  surface  is  proportional  to  the 
velocity.  It  has  also  been  shown  (Tully  and  Preston,  197l)that 
the  main  contribution  to  this  derivative  comes  from  the 
velocity  of  the  diatomic  that  is  crossing  the  seam.  Thus, 
if  the  diatomic  crosses  the  seam  in  an  expaneion  mode  this 
would  correspond  to  going  to  the  upper  surface  so  that  P^ 
should  decrease  and  Pj  should  Increase.  If  the  diatomic  is 
contracting  when  going  through  the  seam,  the  system  ends  up 
on  the  lower  surface  so  that  Py  should  increase  and  P, 
should  decrease.  Considering  the  results  for  Pj  (Fig.  57) 
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Fig.  (561  Plots  of  and  vs.  t for  ths  non-adlabatic 
collision  of  h"^  with  Hj  on  ths  lowest  two 
surfaces.  n=0,  Etot“  -1.9175  e^V.,  2^-1. 45  a^^ 
and  P^  «.S98  m„a  t"1. 

p o 
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one  sees  that  F2  increases  rapidly  to  about  .1  for  t around 
3 T.  From  Fig.  (48)  this  corresponds  to  positive  values  of 
^2^ ' BC  then  begins  to  contract  so  that  is  negative  and 
decreases.  This  is  followed  by  another  increase  in  ^2 
about  .3  for  t around  4 x.  From  Fig.  (48)  ?2^  is  seen  to  be 
positive  during  the  increase  and  zero  at  the  maximum.  As  BC 
contracts.  ?2  decreases  and  attains  a final  value  of  about 
.OS.  The  larger  increase  in  Fj  is  seen  to  correspond  to  a 
greater  expansion  of  BC. 

Two  important  inferences  about  trajectories  that  origi- 
nate in  the  n equal  to  aero  vibrational  level  can  be  made. 
First,  non-adiabatic  effects  should  only  be  Important  as  the 
particles  exit  the  collision  region.  This  is  so  because 
before  the  collision  the  oscillations  of  BC  are  small  and  do 
not  enter  the  seam  and  because  the  couplings  are  small  for 
small  values  of  r.  Second,  the  mechanism  that  determined 
whether  or  not  rearrangement  should  take  place  should  apply 
here  too.  This  is  because  the  mechanism  depended  on  the 
initial  part  of  the  collision  where  non-adiabatic  effects 
are  small. 

The  second  set  of  trajectories  to  be  considered  here 
also  originate  in  the  n equal  aero  vibrational  level.  The 
initial  values  of  and  1.25  Sq  and  .4261 

respectively.  The  results  of  r,  p and  0,  Pg  are  shown  in 
Figs.  (58)  and  (59)  respectively.  One  sees  that  the  final 
values  of  9 are  close  to  r so  that  this  is  an  example  of 


a rearrangement  process.  The  larger  oscillations  of  6 after 
the  collision  indicates  an  increase  in  vibrational  energy. 

The  results  of  2^,  and  zi,  are  plotted  in  Figs. 
(SO)  and  (61)  respectively.  One  notices  that  during  the 
initial  part  of  the  collision  SC  is  expanding,  which  based 
on  previous  discussions  results  in  a rearrangement.  The 
results  for  and  aj,  are  shown  in  Figs.  (62)  and 

(63)  respectively.  As  seen  in  Fig.  (62), becomes  constant 
after  the  collision  and  this  was  used  as  the  criterion  to 
terminate  the  calculation. 

W is  shown  in  Fig.  (66)  and  the  diabatlc  potentials  Ej 
and  Ej  are  shown  in  Figs.  (65)  and  (66)  respectively.  One 
sees  that  there  is  considerable  vibrational  excitation  in  W 
and  e-i  after  tiie  collision.  One  also  notices  that  Ej  and  E2 
are  close  to  each  other  at  around  t equal  to  6 r and  it  will 
be  seen  that  and  Pj  undergo  large  changes  in  this  region. 

The  results  of  and  Pj  are  shown  in  Fig.  (67).  From 
Figs.  (61)  and  (63)  one  sees  that  the  small  initial  increase 
in  ?2  is  due  to  the  expansion  of  ac.  Then  as  A approaches  B 
it  passes  through  the  seam  and  Pj  decreases  to  very  small 
values.  For  t between  4 t and  5 t the  particles  are  close 
together  and  the  couplings  are  small  so  that  Pj  remains  small. 
At  around  S t Pj  begins  to  sharply  increase  and  from  Figs. 

(61)  and  (63)  this  is  seen  to  result  from  both  AB  and  BC 
expanding.  A passes  through  the  seam  and  exits  and  a 
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Pig  (58)  Plots  of  r and  P vs.  t for  the  non-adiahatio 
collision  of  H with  Hj  on  the  lowest  two  H* 
surfaces.  n«0.  -1.9175  e.V.,  *^=1.25  a 

and  -.4261  mpagt-l. 


Fig.  {591  Plots  of  8 and  Pg  vs.  t for  the  non-adiabatic 
collision  of  H with  on  the  lowest  two  Hj 
surfaces.  n»0,  £tot~  “1.9175  e.v.,  sj*1.25  a^ 
and  P‘^-.1261 
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Fig.  (64)  Plot  of  W vs.  t for  the  non-adiabatlc 

collision  of  with  Hj  on  the  lowest  two 
surfaces.  n=0,  -1.9175  e.V.,  »^=1.25  a^ 

and  p‘  = .4261  «pajjt-l. 
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Pig-  (65}  Plot  o£  the  diabetic  potential  vs.  t for 
the  collision  of  H*  with  Hj  on  the  lowest  two 
Hj  surfaces,  n-0,  £^34=  -1.917S  e.V., 

*}»1.2S  a and  P^  =.4261  m a 
o *1  p o 
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Pig.  (66)  Plot  o£  the  diabetic  potential  vs.  t for 
the  non-adiabatic  collision  of  H'*’  with  Hj  on 
the  lowest  two  surfaces,  n*0,  Btot-  “1.9175 

e.V.,  1^=1. 25  a and  pi  -.4261  m a t*!. 

1 Q tj  p o 
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decreases 


coincide  with  a zero  of  ^ contracts  P2 

and  becoiaes  constant  when  AB  is  away  from  the  seam.  There 
is  a final  increase  in  ?2  at  around  8 t and  this  is  due  to 
an  AB  expansion.  The  reason  that  this  is  not  followed  by 
another  decrease  in  P2  is  that  as  r gets  larger  the  seam 
becomes  narrower  so  that  by  the  time  AB  has  fully  expanded 
it  is  no  longer  in  the  seam. 

The  final  two  sets  of  trajectories  to  be  discussed  here 
originate  in  the  n = 4 vibrational  level.  From  Eqs.  (6-33) 
and  (6-34)  one  finds  that  is  equal  to  -1.15  mgSg  . The 

values  of  Sj  and  for  the  first  set  of  trajectories  were 
chosen  to  be  2.17  a^  and  .635  mpS^  respectively. 

The  results  of  r,  P^  and  6,  P^  are  shown  in  Figs.  (68) 
and  (69)  respectively.  As  is  to  be  expected,  the  initial 
oscillations  in  P^.,  6.  and  Pg  are  noticeably  larger  than  in 
the  previous  examples.  One  sees  that  the  final  rate  of 
increase  of  r is  leas  than  its  initial  rate  of  decrease  so 
that  the  internal  energy  of  BC  has  increased.  However, 
the  final  oscillations  of  Pg  are  smaller  so  that  the 
increase  is  probably  due  to  an  electronic  transition. 

The  results  of  2^,  Pj^  and  z^'  ®zi  shown  in  Figs. 
(70)  and  (71)  respectively.  The  invariance  in  the  final 
values  of  Pg^  was  used  to  terminate  the  calculation.  One 
notices  that  as  A and  B begin  to  Interact  at  around  4 z 
BC  is  contracting  which  again  leads  to  non-rearrangement. 
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Fig,  (60)  Plots  of  r and  P vs,  t for  the  non-adlabatic 
collision  of  U with  ^*0  lowest  lit 

surfaces,  n=4,  ^tot”  ' 


.9175 


Fi?. 


Plots  of  9 and  P3  vs.  t for  the  non-adlabatlc: 
collision  of  with  Hj  on  the  lowest  two  a! 
surfaces,  n-9,  E^ot^  -1.9175  e.v. , sJ-2.17  a 
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Ill 

Zj  and  Pjj  are  presented  In  Fig.  (72)  and  will  be  used  later 
to  help  understand  the  behavior  of 

W is  shown  in  Fig.  (73)  and  one  sees  that  after  the 
collision  there  is  a substantial  increase  in  W with  smaller 
oscillations  which  indicates  that  the  motion  is  taking 
place  on  the  upper  surface.  The  diabatic  potentials  and 
£2  shown  in  Figs.  (74)  and  (75)  and  one  notices  that 
there  are  crossings  between  and  at  around  5 t. 

The  results  of  and  P2  are  given  in  Fig.  (76)  . One 
notices  a considerable  increase  in  structure  but  some 
understanding  can  be  achieved  by  considering  the  general 
trends.  From  Fig.  (71),  one  sees  that  the  initial  increase 
in  ?2  i®  ^ an  expansion  of  Be.  ?2  decreases  somewhat 
as  BC  contracts  and  then  begins  oscillating  rapidly  between 
value  of  lero  and  .6.  From  Figs.  (71)  and  (72)  one  sees 
that  the  oscillations  increase  as  BC  expands  across  the 
seam  and  then  decrease  as  A crosses  the  seam.  The  increase 
in  the  next  set  of  oscillations  is  due  to  A crossing  the 
seam  as  it  exits  and  this  followed  by  a decrease  due  to  the 
contraction  of  BC.  As  BC  expands  P2  starts  to  oscillate 
between  zero  and  one  and  attains  a value  of  about  0.6  when 
SC  has  contracted  out  of  the  seam.  As  in  the  previous 
example  there  is  a final  increase  in  Pj  as  BC  again  expands 
into  the  narrowing  seam. 

Thus  even  though  there  is  an  enormous  amount  of 
structure  in  the  results  for  the  probabilities  their 


222 


223 


224 


the  non-adiahatic  collision  of  H with  on 
the  lowest  two  surfaces,  n=4,  -1,9175 

e,V,,  *J-2,17  a^  and  p‘  -.635 
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Plot  of  the  diabetic  potential  E 
non-adiabatic  collision  of  h'*’  m1 
lowest  two  surfaces  of  Hj.  n-4, 
e.V.,  Zj-2.17  a and  =.635  m_ 


Fl9  (76)  Plots  of  pj  and  P^  vs.  t for  the  non-adlabatic 
collision  of  H*  with  Hj  on  the  lowest  two  H* 
surfaces.  n=4,  -1.9175  e.v.,  sf-2.17  a 

and  -.SIS  .n  a --1  ^ ° 


general  behavior  can  still  be  understood  in  terms  of  simple 
seam  crossings.  One  also  notices  that  the  non-adiabatic 
effects  before  the  collision  are  substantially  larger  than 
in  the  previous  n = 0 cases.  This  somewhat  puts  into 
question  whether  the  mechanism  that  determined  whether  or  not 
rearrangement  would  occur  is  valid  here.  However  from  the 
cases  that  have  been  studied,  it  still  seems  to  work. 

The  final  set  of  trajectories  also  start  in  the  n 
equal  4 vibrational  level.  The  initial  values  of  Zj  and 
are  2.17  a,  and  -.435  mpao  reapectively.  The  results  for 

r,  and  8,  are  shown  in  Figs.  (77)  and  (78).  One 
notices  a substantial  increase  In  the  absolute  value  of 
after  the  collision  which  indicates  that  internal  energy 
has  been  lost.  Thus  one  would  expect  that  the  probability 
for  a transition  to  the  upper  surface  to  be  small.  One 
notices  that  the  final  values  of  8 are  close  to  x so  that 
this  is  an  example  of  a rearrangement  process. 

Fig.  (79)  displays  the  results  tor  zj,  and  Pj^.  The 
only  difference  in  boundary  conditions  between  this  example 

and  the  previous  one  is  the  sign  of  Pi  . Based  on  the 
*1 

mechanism  introduced  earlier,  one  would  predict  that  a 
rearrangement  would  occur  and  it  is  gratifying  to  see  that 
it  indeed  does.  The  results  for  z^  and  P^  are  shown  in 
Pig.  (80)  and  will  be  used  later  when  the  probabilities  are 


discussed. 


Fig.  (771  Plots  of  r and  vs.  t for  the  non-adiabatic 
collision  of  H*  with  Hj  on  the  lowest  two  B* 
surfaces.  n=4,  -1.9175  e.V.,  zj=2.17  a 


Fig. 


(78)  Plots  of  9 and  Pg  vs.  t for  the  non-adiabatic 
collision  of  H with  on  the  lowest  two 
surfaces.  n=4,  Etot=  -1.917S  e.V.,  s^-2.17  a^ 
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Fig.  <79)  Plots  Of  and  vs.  t for  the  non-adiabatlc 
collision  of  with  on  the  lowest  two  H* 
surfaces.  n=4,  Etot"  -1.9175  e.V. , *J=2.17  a 
and  P^  - -.635  m a^T'l.  ° 
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Fig.  (81) 


collision  of  H with 
suxfaces.  n*4,  E. 


non- adiabatic 
on  the  lowest  two  nt 


-1.9175 
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W is  shown  in  Pig.  (81)  and  one  sees  that  after  the 
collision  W undergoes  small  oscillations.  Thus  one  would 
expect  a small  final  value  of  Pj. 

and  Pj  are  shown  in  Fig.  (82)  and  again  there  is 
considerable  structure,  From  Fig.  (79)  one  sees  that  the 
initial  increase  in  Pj  expansion  of  BC.  BC 

then  contracts  and  P^  remains  constant.  As  BC  again  expands 
across  the  seam  Pj  starts  to  oscillate  between  small  values 
and  values  that  are  close  to  0.7.  The  oscillations  in  P2 
begin  Co  decrease  as  A crosses  the  seam.  The  second  set  of 
oscillations  are  due  to  an  expansion  of  AB  followed  by  a 
contraction  of  BC.  The  oscillations  are  small  since  the 
couplings  in  this  region  are  small.  The  increase  in  the 
final  set  of  oscillations  is  due  to  C crossing  the  seam  as 
It  exits.  The  oscillations  decrease  as  AB  contracts  and 
Pj  attains  a final  value  of  about  0.1. 

As  a final  comment  one  notices  that  the  non-adiabatic 
effects  are  much  more  pronounced  in  the  Initial  part  of 
the  collision.  This  gives  support  to  the  mechanism  that 
predicts  rearrangement  since  its  prediction  is  correct  even 
when  the  non-adiabatic  effects  are  large. 

6-3  Experimental  Studies 

The  earlier  experimental  studies  involving  the  colli- 
sion of  a proton  with  H2  focused  on  obtaining  probabilities 
of  electron  transfer  as  a function  of  the  Initial  energy 
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of  the  proton  at  a fixed  laboratory  angle  IZlemba,  Lockwood, 
Piorgan  and  Everhart,  1960)  . The  range  of  energy  of  the 
proton  was  between  one  and  200  XeV  and  the  laboratory 
angle  was  in  the  neighborhood  of  5* . Por  such  large  pro- 
ton energies  this  angle  corresponds  to  almost  head  on 
oollisions.  The  study  done  by  Ziecnba  et  al.  (Ziemba, 
Loclcwood,  Morgan  and  Everhart,  1960}  showed  that  the  system 
exhibited  a resonant  electron  capture,  i.e.  the  probability 
oscillated  around  a value  of  about  .4.  When  the  probability 
was  plotted  versus  the  reciprocal  velocity  it  was  found  that 
the  spacings  between  adjacent  resonances  were  almost  the 
same.  A comparison  of  the  probability  of  electron  transfer 
for  the  systems  h'*',  and  H , H was  done  by  Loc)cwood  and 
Everhart  (Locliwood  and  Everhart,  19621  and  their  results 
showed  that  the  resonances  for  the  two  systems  occurred  at 
nearly  the  same  energies.  The  positions  of  the  resonances 
for  the  system  H , H is  well  understood  and  can  be  simply 
related  to  the  nximber  of  times  the  electron  jumps  freon  one 
proton  to  the  other.  This  suggests  a similar  mechanism  for 


The  next  series  of  experiments  to  be  considered  here 
used  a much  lower  kinetic  energy  for  the  incoming  ion.  The 
study  by  Krenos  and  Wolfgang  (Krenos  and  Wolfgang,  1970) 
measured  the  relative  yields  in  the  forward  direction  of 
the  reactive  products  {products  corresponding  to 
raent,  electron  transfer  and  electron  transfer  wit 


rearrange- 


cearrangenent)  with  the  initial  kinetic  energy  of  the  ion 
between  X and  B e.V.  As  expected/  their  results  showed  that 
below  the  threshold  for  molecular  ion  formation,  about 
l.S  e.V.,  reactive  products  due  to  rearrangement  only  were 
observed  and  their  relative  yields  decreased  as  the  energy 
increased  up  to  about  6 e.V.  As  the  energy  increased  above 
threshold,  the  relative  yields  for  electron  transfer  proces- 
ses increased  sharply  and  then  began  to  decrease  in  the 
neighborhood  of  3 e.V,  Between  2 and  5 e.V.  the  relative 
yields  for  the  electron  transfer  processes  were  comparable 
to  the  one  for  rearrangement  only  and  in  some  cases  were 
even  dominant. 

One  of  Che  hopes  of  this  experiment  was  to  find  a simple 
mechanism,  i.e.  complex  formation  or  direct  interaction, 
that  would  acoount  for  the  results.  In  the  study  of  D* 
colliding  with  HD  it  was  found  that  the  ratio  of  the 
relative  yield  of  HD"^  to  Dj*  was  roughly  2:1.  This  is  what 
one  would  expect  if  a long  lived  complex  were  formed. 

However  in  the  study  of  colliding  with  D2  the  relative 
yield  of  D2*  was  about  twice  that  of  HD*  which  indicated  a 
more  direct  mechanism.  This  led  to  their  conclusion  that 
only  one  simple  mechanism  would  not  be  able  to  account  for 
all  of  the  results.  Another  observation  worth  noting  is 
that  there  was  a large  tendency  to  transfer  energy  from 
to  vibrational  modes . 
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similar  results  were  obtained  by  Holliday  et.  al. 
(Holliday,  Huokerman  and  Friedman,  1971).  They  measured 
cross  sections  for  the  reactive  products  of  collisions  of 
D*  with  Hj,  H'*  with  Dj  and  □’’with  HD  with  the  incident  ion 
energy  between  1 and  0 e.V.  They  again  found  that  some 
results  indicated  the  formation  of  a complex  while  others 
seemed  to  involve  direct  processes.  Their  conclusions  were 
that  complex  formation  was  important  for  incoming  ion  ener- 
gies less  that  4 e.V,  but  that  direct  processes  dominated 
at  higher  energies.  They  alec  found  a large  tendency  to 
transfer  energy  from  translational  to  vibrational  inodes. 

Cross  sections  for  the  reactive  products  resulting  from 
the  collisions  between  H*  and  Dj  with  the  incident  proton 
energy  between  .3  and  100  e.V.  were  measured  by  Haier 
(Maier,  1971).  The  Increased  range  of  proton  energy  allows 
the  investigation  of  the  effect  of  the  breakup  channel  which 
has  a center  of  mass  threshold  of  about  5 e.V.  They  found 
that  the  total  reaction  cross  section  was  much  smaller  tlian 
the  Langevin  cross  section  (Gioumousis  and  Stevenson,  1958). 
The  values  smaller  that  theoretically  predicted  were  also 
observed  by  the  previous  authors. 

Also  obtained  in  this  study  were  the  branching  ratios 
tor  the  reactive  products.  The  branching  ratio  is  defined  as 
the  cross  section  for  the  particular  reactive  product  divided 
by  the  sum  of  the  cross  sections  for  all  reactive  products. 

As  was  expected,  the  branching  ratio  for  D below 
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the  threshold  for  molecular  ion  formation  was  one.  Above 
threshold  it  sharply  decreased  and  reached  a minimum  at  a 
center  of  mass  energy  of  about  S e.v.  It  then  increased  by 
a factor  of  three  between  5 and  10  e.V.  and  remained  fairly 
constant  above  10  e.v.  This  increase  is  probably  due  to 
formed  through  a breahup  process.  The  branching  ratio  for 
□2  is  zero  for  center  of  mass  proton  energies  less  than 
about  2 e.v.  It  increases  rapidly  between  2 and  S e.v.  and 
levels  off  at  a value  greater  than  .5  for  energies  greater 
than  5 e.V.  Thus  electron  transfer  is  the  dominant  process 
above  5 e.v.  The  branching  ratio  for  HD*  rises  sharply  as 
the  energy  is  increased  above  threshold.  It  reaches  a 
maximum  in  the  neighborhood  of  2.5  e.v.  and  guichly  decreases 
to  small  values  as  the  energy  Increases.  This  is  further 
evidence  Chat  complex  formation  is  important  for  low  energies 
but  that  at  higher  energies  direct  processes  dominate. 

A rather  extensive  study  of  this  system  with  the 
collision  energy  between  1 and  8 e.V.  was  done  by  Krenos 
et  al.  (Krenos,  Preston,  Wolfgang  and  Tally,  1974) . They 
studied  the  collisions  of  H*  with  D2,  D*  with  HD  and  D* 
with  D2  and  obtained  the  cross  sections  for  the  various 
reactive  products.  They  compared  these  with  those  obtained 
in  the  previous  studies  and  found  that  there  was  fairly 
good  agreement. 

A crossed  molecular  beam  arrangement  was  used  in  this 
study  to  obtain  a velocity  analysis  of  the  products  at 
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several  lab  angles.  These  results  were  compared  Co 
theoretical  calculations  (to  be  discussed  in  the  next 
section)  to  give  a better  understanding  of  the  mechanics 
that  govern  the  reactions  as  a function  of  the  initial 
collision  energy.  In  general  they  found  that/  for  the 
initial  diatonic  in  the  lowest  vibrational  level,  non- 
adiabatic  effects  were  not  important  as  the  particles 
approached  each  ocher.  This  is  because  the  vibrations  of 
the  diatomic  ace  not  large  enough  to  reach  the  avoided 
crossing.  As  the  particles  get  close  together  they  enter 
the  deep  well  region  of  the  lower  surface.  Since  the 
surfaces  are  far  apart  in  this  region  non-adiabatic  effects 
are  again  unimportant  so  that  the  nuclear  motions  take 
place  on  the  lowest  surface.  This  well  region  was  respon- 
sible for  complex  formation  and  for  the  conversion  of 
kinetic  to  internal  energy.  As  the  particles  separated, 
non-adiabatic  effects  became  important  and  in  some  cases 
there  were  numerous  crossings  of  the  seam. 

Some  insights  into  the  mechanieras  that  led  to  final 
nuclear  configurations  were  also  obtained.  It  was  found 
that  above  a center  of  mass  collision  energy  of  4.S  e.V, 
complexes  no  longer  formed.  This  was  expected  since  the 
dissociation  energy  of  the  initial  diatcmic  is  4.5  e.V.  and 
for  a collision  energy  above  this  there  is  no  way  that 
energy  can  be  shared  among  the  three  particles  to  give  a 
stable  complex.  It  was  shown  that  above  7 e.V.  velocity 
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contours  fron  theoretical  calculations  using  hard  sphere 
potentials  agreed  fairly  well  with  experimental  ones.  This 
indicated  that,  in  this  region  of  energy,  nuclear  motion  is 
governed  mostly  by  repulsive  forces  and  that  direct  processes, 
i.e.  knockout,  etc.,  led  to  rearrangement.  The  theoretical 
calculations  also  showed  that  as  the  energy  decreased  below 
4.5  e.V.  Che  collision  time  increased  which  indicated  that 
complex  formation  was  important. 

A further  experiment  done  by  Ochs  and  Teloy  (Ochs  and 
Teloy,  1974)  focused  on  obtaining  more  accurate  cross  sections 
for  the  reactive  products.  The  interesting  aspect  of  this 
work  was  that  the  agreement  with  previous  experimental 
studies  was  only  fair  but  the  agreemsnt  with  the  theore- 
tical results  was  excellent. 

The  crossed  beam  experiment  done  by  Lees  and  Rol  (Lees 
and  Rol,  1975)  concentrated  on  obtaining  energy  distributions 
for  the  reactive  products.  Their  agreement  with  the  theo- 
retical predictions  was  also  excellent.  The  added  feature 
in  this  experiment  was  that  the  initial  diatomic  could  have 
a vibrational  quantum  number  of  up  to  seven.  Their  results 
indicated  that  the  higher  vibrational  levels  noticeably 
enhanced  the  formation  of  molecular  ions. 

6-4  Theoretical  Studies 

Since  the  system  has  only  two  electrons,  it  is 
possibls  to  obtain  the  potential  energy  surfacds  by  using 
ab  initio  methods.  This  is  of  particular  importance  in 
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trajectory  calculations  because  comparing  experimental 
results  to  ab  initio  theoretical  results  should  give  some 
indication  of  the  validity  of  using  classical  trajectories 
as  a tool  to  study  chemical  reactions.  This  was  one  of  the 
primary  incentives  for  the  calculations  done  by  Csizmadia 
et  al.  (csizmadia,  Polanyi,  Roach  and  Wong,  1969),  They 
calculated  differential  cross  sections  for  the  production  of 
H*  from  the  reaction  of  D'*'  with  H2  Initial  center  of  mass 
relative  energies  of  1 and  4.5  e.V.  The  three  dimensional 
classical  trajectories  were  calculated  using  the  lowest 
adiabatic  surface  which  was  determined  by  ab  initio  means. 
When  compared  to  experimental  results  it  was  found  that  the 
cross  sections  were  too  large.  This  was  not  too  surprising 
since  this  calculation  neglected  the  energetically  open 
channels  corresponding  to  electron  transfer  processes. 

This  wort  was  important  however  because  it  gave  motivation 
for  the  experiments  to  be  done. 

A calculation  referred  to  in  the  previous  section  that 
incorporated  the  upper  surface  was  done  by  Tully  and 
Preston  (Tully  and  Preston,  1971) . Comparison  to  experimental 
results  were  also  made  (Krenos,  Preston,  Wolfgang  and  Tully, 
1974)  . The  theoretical  approach,  termed  the  Trajectory 
Surface  Hopping  model,  handled  the  dynamics  of  the  reaction 
in  the  following  way.  Initially,  the  diatomic  (D2)  was  in 
a low  vibrational  level  so  that  seam  crossings  did  not  occur 


prior  to  the  collision.  In  this  region,  nuclear  motion 
evolved  on  the  lower  surface  only.  When  the  incoming  ion 
(H'*’)  closely  approaches  the  diatomic,  the  particles  are  in 
the  deep  well  region  of  the  lower  potential  and,  since  the 
couplings  are  small,  nuclear  motion  is  again  governed  by 
the  lowest  surface.  As  the  particles  separate,  seam  cros- 
sings can  occur  so  that  non-adiabatic  effects  must  be 
included.  If  a seam  crossing  occurs  the  trajectory  splits 
into  two  trajectories  corresponding  to  motion  on  each 
surface.  The  trajectories  are  assigned  probabilities  that 
are  determined  from  the  non-adiabatic  couplings.  These 
couplings  are  rather  easy  to  obtain  since  the  surfaces  are 
obtained  from  a Oiatomics  in  Molecules  approach.  Further 
seam  crossings  result  in  further  splitting  of  the  trajec- 
tories so  that  a unique  initial  trajectory  usually  evolves 
into  a number  of  weighted  final  trajectories  all  of  which 
have  unique  final  conditions.  Cross  sections  are  obtained 
by  a Monte  Carlo  averaging  over  the  initial  conditions. 
Other  details  such  as  the  assignment  of  positions  and 
momenta  after  a surface  crossing  has  occurred  are  discussed 
in  the  articles  referred  to  above. 

The  excellent  agreement  between  theory  and  experiment 
was  very  encouraging  since  it  showed  that  trajectory  calcu- 
lations could  still  be  used  for  multi-surface  systems  if 
the  couplings  could  be  incorporated  in  a reasonable  way. 
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As  is  apparent  from  the  previous  sectionr  this  is  important 
because  trajectory  calculations  combined  with  experimental 
data  often  lead  to  a better  understanding  of  the  dynamics  of 
the  reaction. 

A quantum  mechanical  calculation  for  the  collinear  Hj* 
system  was  done  by  Top  and  Baer  (Top  and  Baer,  1977).  They 
used  a diabatic  electronic  representation  and  the  total 
energy  was  in  the  few  e.v.  range.  Their  study  showed  that 
for  a given  total  energy,  the  probability  for  electron 
transfer  was  greatly  enhanced  if  the  diatomic  was  initially 
in  a higher  (^4)  vibrational  level.  Since  the  n « 4 
vibrational  level  of  Hj  is  nearly  coincident  in  total 
internal  energy  with  the  n = 0 vibrational  level  of  H2*' 
this  behavior  was  attributed  to  a resonance  between  these 
two  levels.  For  lower  initial  vibrational  levels  they 
argued  that  the  deep  well  in  the  lower  surface  caused 
vibrational  excitation  which  populated  the  upper  vibrational 
levels.  As  the  particles  leave  the  reaction  zone  resonances 
between  these  populated  levels  and  those  of  the  H2*  lead  to 
electronic  transitions.  We  shall  return  to  some  of  the 
more  numerical  aspects  of  this  study  in  the  next  section. 
d-5  Comparisons 

Since  this  work  only  considers  the  collinear  confi- 
gurations of  the  Hj  system  no  direct  quantitative  comparisons 
can  be  made  with  either  the  Trajectory  Surface  Bopping 
studies  or  the  experimental  studies  discussed  in  Section  6-3. 
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It  is  encouraging  to  see  however  that  some  of  their  conclusions 
are  borne  out  in  this  study  also.  The  oalculations  discussed 
in  section  6-2  also  indicate  that  non-adiabatic  effects  are 
most  important  as  the  particles  are  leaving  the  collision 
region.  In  other  words,  the  nuclear  motions  before  the 
oollision  and  while  the  particles  are  all  close  together  are 
governed  by  the  lowest  surface.  The  simple  mechanism  sug- 
gested in  section  fi-2a,  that  determines  whether  rearrangements 
take  place  for  the  adiabatic  test  case,  seems  to  hold  for  the 
two-surface  case  also.  This  supports  the  conclusions  dis- 
cussed in  sections  6-3  and  6-4  that  the  lowest  surface 
determines  whether  rearrangements  occur.  In  concurrence  with 
the  experimental  and  trajectory  studies,  the  results  of 
Section  6-2  also  indicate  that  there  is  a large  tendency  to 
transfer  energy  between  vibrational  and  translational  modes. 
Finally  t-he  results  of  Section  6-2  indicate  that  they  can 
not  be  toatally  explained  in  terms  of  either  direct  mech- 
anisms or  complex  formation.  This  is  not  unexpected  since 
the  total  energy  used  in  this  study  is  between  the  higher 
energies  where  direct  processes  dominate  and  lower  energies 
where  complex  formation  is  important. 

The  quantum  mechanical  calculation  by  Top  and  Baer 
determined  the  total  electronic  transition  probabilities  for 
a number  of  total  energies.  As  discussed  in  Section  3-7, 
it  is  possible  to  construct  total  electron  transfer  proba- 
bilities using  this  formalism  if  appropriate  averages  are 
tahen  over  the  initial  conditions. 


For  the  systein,  the  initial  conditions  for  the 
diatonic  that  need  to  be  averaged  over  are  the  vibrational 
coordinate  and  momenta  of  the  initial  ^2  diat<»iic.  The 
vibrational  distributions  for  the  coordinate  are  propor- 
tional to  the  eigenfunctions  of  the  Horse  potential  (Horser 
1929) . The  season  why  they  are  only  proportional  to  these 
functions  is  that  some  of  the  quantum  mechanical  distribution 
lies  outside  of  the  classical  turning  points  which  is 
physically  meaningless  in  this  formalism.  Instead  one  must 
require  that  the  Horse  eigenfunctions  are  normaiized  between 
the  classical  turning  points.  This  is  easily  done  by 
determining  the  integral  of  IVin!*  between  the  turning  points, 
where  is  the  Morse  eigenfunction  for  the  n**  vibrational 
level.  The  normalization  constant  is  then  defined  as  the 
inverse  square  root  of  this  integral. 

There  still  remains  the  problem  of  averaging  over  the 
initial  momentum  of  the  diatceiic.  In  actuality  the  charge 
transfer  probability  should  be  found  by  integrating 
*’2*‘l^®''®n I *n***  I * I ’ over  the  curve  of  constant 

in  classical  phase  space.  However  it  was  assumed  for 
simplicity  that,  for  a given  value  of  the  vibrational 
coordinate,  tne  corresponding  momentum  may  be  either  positive 
or  negative  with  equal  probability.  With  this  assumption, 
the  total  electron  transfer  probability  is  defined  as 


;d*p-.^(E,E„,z.a)|^'„(z)|2, 
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Each  equally  spaced  value  of  z corresponda  to  an  integer 
value  of  y between  1 and  N . 

The  results  for  the  integrands  and  in 
Eq.  (6-44)  for  n = 0 to  n = 4 are  shown  in  Pigs.  (83)  through 
(87)  respectively.  In  each  case  50  equally  spaced  values  of 
z were  used  so  that  y goes  fron  1 to  50.  A logarithmic 
scale  for  the  probabilities  was  used  in  order  to  show  more 
of  their  behavior  when  the  probabilities  are  small.  The 
points  labeled  with  a a correspond  to  cases  where  dissociation 
of  the  complex  did  not  occur  within  a time  of  50  t.  When 
necessary,  lines  between  adjacent  points  were  drawn  so  that 
the  behavior  of  the  probabilities  would  be  easier  to  follow. 

All  calculations  were  done  at  a total  energy  of  -1.918  e.V. 
and  the  initial  conditions  were  determined  as  in  Section  6-2. 
The  R|N  notation  significies  that  one  is  going  from  a region 
where  rearrangements  occur  to  a region  where  non-rearrangements 
occur.  N|r  has  a similar  meaning. 

One  sees  in  Figs.  (83)  through  (35)  that,  for  n<2,  there 
are  two  fairly  localized  regions  where  the  proloability  is 
larger  than  0.1.  As  is  seen  in  Figs.  (86)  and  (87)  this 
is  no  longer  the  case  for  n>3.  One  notices  that  the  regions 
where  the  probability  is  greater  than  0.1  are  considerably 
more  extensive  for  n^3  which  should  give  rise  to  substan- 
tially larger  total  electron  transfer  probabilities.  A 
rough  estimate  of  the  behavior  of  the  total  electron  transfer 
probability  as  a function  of  n can  be  made 'by  simply  counting 


the  number  of  times  the  integrand  in  Eq.  (6-44)  is  greater 

greater  than  0.1  for  the  Integrand  23,  16,  16,  32  and  57 
times  respectively.  This  indicates  that  the  total  electron 
transfer  probability  should  not  change  very  much  for  n<2 
and  that  it  should  noticeably  increase  for  n>3. 

Total  electron  transfer  probabilities  were  obtained  by 
integrating  the  results  shown  in  Figs.  (83)  through  (87) 
using  the  trapezoidal  rule.  The  final  probabilities  for 
trajectories  that  did  not  dissociate  within  50  t,  i.e. 
points  labeled  by  d in  the  figures,  were  assigned  a value  of 
0,5  + 0.5.  The  results  of  this  work  are  compared  to  those 
estimated  from  Pig.  (11)  in  the  paper  by  Top  and  Baer  (Top 

trajectories  that  did  not  dissociate.  The  agreement  for 
n * 0 and  4 is  quite  good  and  well  within  the  estimated 
errors.  For  n ■ 1 the  agreement  is  not  as  good  but  the 
results  are  still  within  40%  of  each  other.  It  is  encour- 
aging to  see  that  both  sets  of  results  show  relatively 
little  change  in  the  probability  for  n going  from  0 to  1. 

It  ie  also  pleasing  to  see  that  both  studies  show  a sharp 
increase  in  the  probability  for  n = 4.  Sources  of  disagree- 
ment could  arise  from  excluding  the  region  of  configuration 
space  beyond  the  turning  points  and  from  not  using  enough 
initial  values  of  the  diatomic  coordinate.  The  approxi- 
mations used  in  obtaining  Eq.  (6-43)  could  also  be  a source 
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b)  Estimated  from  F19,  (11)  in 
[Top'and  Baer,  1977) 


of  error.  Ml  of  these  sources  of  error  cen  be  numerically 
checked  and  could  be  the  subject  of  future  studies. 

The  most  noticeable  disagreement  in  the  results  occurs 
for  n a 3.  This  was  called  an  "intermediate  case"  (inter- 
mediate in  the  sense  that  the  mechanism  for  charge  transfer 
is  neither  totally  due  to  resonances  for  n>4  nor  to 
E>opulating  the  upper  vibrational  levels  via  the  deep  well 
for  n<2)  in  the  quantum  mechanical  study  and  was  not  further 
considered.  The  reason  for  such  a low  probability  is  not 
clear.  It  would  seem  that  an  "intermediate  case"  would  lead 
to  results  somewhere  between  the  ones  for  n<2  and  the  ones 
for  n>4. 

6-6  Conclusions 

This  work  primarily  addressed  two  fundamental  problems 
in  scattering  theory.  One  of  the  problems  centered  around 
the  electronic  basis  used  in  the  expansion  for  tlie  full 
solution.  The  other  problem  focused  on  how  to  handle  the 
problem  for  the  nuclear  degrees  of  freedom  when  mere  than 
one  electronic  basis  function  is  included  in  the  expansion. 

In  chapter  two  it  was  shown  that  the  Hamiltonian  for 
the  nuclear  degrees  of  freedom  was  not  hermitian  in  the 
adiabatic  electronic  representation.  This  is  of  little 
consequence  when  only  one  electronic  state  is  energetically 
accessible.  However,  when  it  oecame  necessary  to  include 
mare  than  one  electronic  state  it  was  pointed  out  that 
other  electronic  representations  could  be  more  suitable. 
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The  diebatic  and  alaosc  adiabatic  representations  discussed 
in  sub-sections  2-2b  and  2-2c  respectively  were  shovm  to 
restore  hermiticity  to  the  Hamiltonian.  The  minimisation 
procedure  discussed  in  section  2-3  was  shown  to  lend  to 
well  behaved  potentials  but  its  usefulness  depended  on 
whether  the  coupling  terms  after  transformation  could  be 
neglected  or  not. 

Given  a suitable  electronic  basis«  the  nuclear  problem 
was  the  subject  matter  of  chapter  three.  It  was  shown  that 
if  the  nuclear  expansion  coefficients  for  the  time  independ- 
ent SchrSdinger  eqxiation  are  vritten  in  the  form  of  a common 
eikonal  and  the  short  wavelength  approximation  is  used,  a 
transformation  from  coordinate  to  the  time  variable  led  to 
a set  of  coupled  first  order  differential  equations  that 
determined  the  nuclear  positions,  momenta  and  expansion 
coefficients.  It  was  seen  that  this  set  of  coupled  equations 
is  formally  equivalent  to  Hamilton's  equations  of  motion 
where  the  independent  variables  now  Include  the  real  and 
imaginary  parts  of  the  nuclear  expansion  coefficients. 

One  of  the  principal  advantages  that  this  formalism  provided 
was  that  lengthy  expansions  in  internal  nuclear  states  were 
avoided. 

An  application  to  the  oollinear  H3*  system  was  made 
in  chapters  four  through  six.  The  trajectories  presented  in 
section  6-2  indicated  how  useful  they  can  be  in  obtaining 
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a better  understanding  o£  the  basic  inechanisms  Involved  in 
chemical  reactions.  It  was  very  encouraging  to  see  that 
many  of  the  conclusions  about  the  mechanisms  of  the  H***  + 
reaction  made  from  the  experimental  and  Trajectory 
Surface  Hopping  studies  were  supported  by  this  work  also. 

The  reasonable  agreement  with  the  quantum  mechanical  study 
indicated  that  this  approach  is  capable  of  yielding  physi- 
cally meaningful  results  even  when  light  nuclear  masses  are 
involved. 

hs  was  pointed  cut  in  chapter  three,  one  of  the  conse- 
quences of  using  a common  eikonal  is  that  the  final  trajec- 
tories evolve  on  an  average  potential  energy  surface. 
Although  expressions  Such  as  Bn.  (3-108)  allow  one  to 
determine  total  electron  probabilities,  one  still  needs  to 
determine  cross  sections  in  order  to  make  direct  cwiiparisons 
with  experiment.  The  problem  that  arises  is  that  cross 
sections  are  defined  with  quantum  mechanical  initial  and 
final  states  that  evolve  on  their  appropriate  potential 
energy  surfaces.  The  corresponding  trajectories  should 
evolve  initially  and  finally  on  only  one  surface.  As  was 
seen  in  Section  6-2  this  was  rarely  the  case, 

This  problem  could  in  principle  be  solved  by  finding  an 
expression  that  relates  the  initial  conditions  to  the  final 
conditions.  The  initial  and  final  conditions  correspond  to 
points  in  phase  space  that  lie  on  constant  energy  curves 
of  the  initial  and  final  internal  energies  respectively. 
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Once  these  expressions  have  been  found  the  c 
can  be  obtained  In  a straightforward  manner  (Hicha,  1982)  by 
either  conservation  of  flux  or  from  transition  matrix 
elements.  The  cross  sections  obtained  from  transition 
matrix  elements  should  be  more  reliable  since  they  include 
interference  terms. 

In  closing,  it  should  be  pointed  out  that  this  work  is 
rather  preliminary  and  much  more  remains  to  be  done.  In 
particular,  an  extension  to  three  dimensions  would  be 
important  because  it  would  allow  more  direct  comparisons 
to  experiments  and  Trajectory  Surface  Hopping  studies. 

This  would  also  allow  one  to  address  the  problem  of  obtain- 
ing cross  sections.  It  has  also  recently  been  shown  (Micha, 
1962)  that  the  equations  obtained  in  chapter  three  correspond 
to  the  first  term  in  the  expansion  of  the  nuclear  wave 
function  in  powers  of  the  local  deBroglie  wavelength. 
Inclusion  of  the  next  term  in  this  collinear  case  should 
give  some  indication  of  the  error  introduced  by  using  the 
short  wavelength  approximation.  Finally  a 
systems  with  more  than  two  surfaces  would  also  b 
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where  x^,  and  r are  vectors  Crom  the  center  of  mass 
of  the  three  particles  to  particles  Ar  3 and  C respectively 
and  the  channel  subscript  1 has  been  suppressed.  The  tine 
independent  Schrodinger  equation  in  this  coordinate  system 
is  given  by 

5E ’s  ■ 2i  ' 0 (Al-3) 
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and  for  simplicity 
state  is  needed  in 


we  have  assumed  only  one  electronic 
the  6orn*0ppenheimer  expansion.  Letting 
iZfilz)  and  expanding  the  solution  as 
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If  W is  zero  the  solutions  of  Eg.  {Al-E]  could 

be  written  in  terns  of  a product  of  functions  of  z and 
z but  their  eigenvalue  spectrum  would  both  be  continuoiis. 


Defining 
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(sl-8)  becomes 
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end  c and  a are  the  hyperspherical  coordinates.  The 
voluam  element  is 
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Letting 
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One  sees  in  Eg.  (Al-19)  that  R satisfies  an  eqaution 

that  is  formally  equivalent  to  the  two  body  probleu. 

The  essential  difference  is  that  the  potential  in  the 


body  problem  is  considerably 


complicated. 


APPENDIX  2 
COMPUTER  PROGRAM 
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Call 
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Datal  consists  of  tbe  parameters  for  the  model  coupling, 
the  parameters  for  the  diatomic  potentials,  certain 
physical  constants  such  as  the  masses  of  the  particles, 
the  value  of  h,  etc.  and  the  parameters  used  to  terminate 
the  calculation. 

SCCST  is  a subroutine  that  calculates  certain  mass 
dependent  constants  used  throughout  the  program. 

SCTSJ  is  a subroutine  that  determines  the  initial  conditions 
and  calls  the  integration  subroutine. 
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Subroutine  S 


2es 

Data2  contains  Infonnation  for  taking  numerical  deriva- 
tiveSf  parameters  needed  in  the  integration  subroutine, 

DE,  and  information  about  the  initial  state  such  as  its 
total  energy,  the  vibrational  state,  the  initial  relative 
distance,  the  number  of  initial  states,  etc. 

Subroutine  SCINI  determines  the  initial  conditions  for 
the  nuclear  variables  in  terms  of  both  hyperspherical 
and  Jacobi  coordinates.  The  initial  conditions  are 
stored  in  an  array  and  used  in  the  Do  loop. 

Data!  contains  two  integers  K and  H.  If  M«2  and  K^l, 
only  trajectories  with  positive  initial  values  of  the 
diatomic  monentum  will  be  run.  If  M-2  and  K=2  only 
negative  initial  values  will  be  cun.  For  K»1  and  H-1,  all 
trajectories  will  be  run. 

Data4  consists  of  two  sets  of  data.  The  first  set  contains 
arguments  for  time  reversal  and  arguments  used  in  the 
integration  subroutine,  DE,  to  print  intermediate  values 
of  the  trajectories  and  energies.  The  second  set  deter- 
mines which  trajectories  will  be  plotted  In  subroutine 
SCPLOT. 

DataS  are  the  initial  conditions  in  both  hyperspherical 
and  Jacobi  coordinates. 

Subroutine  DE  is  a well  documented  subroutine  [Shampine 
and  Gordon,  1975]  that  integrates  the  equations.  Some 
modifications  were  made  in  order  to  print  intermediate 
the  trajectories  and  energies  and  form  arrays 


values  of 


used  for  plotting.  The  procedure  for  terminating  the 
calculations  discussed  in  Section  6-2  was  also  implemented. 
Subroutine  SCPLOT  plots,  based  on  the  Information  in 
Data4,  any  or  all  of  the  trajectories. 

Data?  contains  the  final  values  of  the  trajectories  for 
the  time  reversed  studies. 
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since  the  integral  over  X contains  all  possible  values  of 
Changing  s would  change  the  values  of  the  integrand 
but  leave  the  integral  unchanged.  To  see  this  assurae  that 


where  the  arguments  and  have  been  suppressed.  The 
point  (X.s.l  can  be  expressed  in  terms  of  by  setting 

X’  - X - a (A3-11) 


■ir* 


Thus  Eq.  {A3-10)  is  satisfied  so  that  the  integral  in 
Eq.  (A3-10)  does  not  depend  on  the  initial  value  of  2. 
Since  the  vibrational  wave  functions  are  normalized.  Eg. 
(A3-4)  becomes 


where  is  the  initial  diatomic  separation.  Instead  of 
averaging  over  X,  one  could  just  as  well  average  over  the 
initial  coordinate  of  the  diatomic  which  gives 


arbitrary  large  initial  relative  distance. 
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